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te INTRODUCTION 


A. OVERVIEW OF THESIS 

This thesis will concentrate on three areas of radar 
Signal processing (RSP) and how they can be simulated to aid 
Student understanding. The primary research done was to 
develop algorithms for the display and manipulation of 
several radar signal processes in order to provide students 
taking radar courses a fuller understanding of the processes 
involved. Aspects of radar signal processing that were 
investigated were the development of the DFT as a filter 
bank for radar situations, the use of windows to create 
weighted DFT filter banks and the generation of pulse 
compression coding schemes under doppler conditions. 

1. The DFT as a Filter Bank 

The background as to how the DFT is used to forma 

filter bank from an input sinusoid is examined. The input 
Sinusoids will be either a complex sinusoid from quadrature 
sampling or a real Sinusoid from single channel sampling. 
The reason only Sinusoids are treated is because it 1s a 
Sinewave carrier that needs to be detected and processed in 
most radar Situations. The student will be able to select 
the pulse repetition frequency (PRF) at which the signal is 


to be sampled and the number of DFT filters that are to be 


considered (N) besides various other parameters. Outputs 
will include a graphical display of a Single DFT Eudes 
response for a range of sinusoidal inputs, an overall DFT 
filter bank response to a range of sinusoidal inputs and the 
Output of the DFT filter bank if a single target were to 
appear at a specified doppler frequency. 
2. Weighted Doppler Filter Bank 

In order to reduce the effects of the sidelobes of 
the filters, the input sequence to the DFT is weighted ina 
various ways that are called windows. The student will be 
able to select the windowing function he wishes to view for 
various input lengths and PRFs. The outputs will be 
displayed graphically versus the output of a rectangular 
window for comparison purposes. 

3. Simulation of Pulse Compression Waveforms 

There will be a brief discussion on the various 
pulse compression modulation schemes available and their 
characteristics. Then the simulation program will be used 
to demonstrate how some codes are doppler tolerant. Codes 
to be generated are the Barker, combined Barker, pseudo 
random and two of the various polyphase schemes (the P3 code 
and the Frank). The student will be able to select the code 
he wishes to observe, the conditions under which he wishes 
to observe it and what doppler effects (if any) are to be 


allowed. Outputs will include a graphical display of the 


compressed waveform and the main parameters associated with 
the pulse compression. In addition a subprogram to show the 
effects of eclipsing loss is available. 

4. Reasons For Simulation 

Older analog radars usually involved no more than a 
transmitter and a receiver connected to a video display. 
This simple and direct system was easy to teach as signal 
processing was at a minimum. However, modern radars are 
much more complex and employ sophisticated radar signal 
processing techniques to extract weak target signals from 
large interference due to rain or ground clutter against 
which the older radars could not discriminate. This type of 
signal processing is very computational intensive and is 
performed in dedicated computers housed within the radar 
system itself. These computers take the sampled signal as 
input and output the finished results such as target 
parameters. Due to the computer’s hardware configuration 
for high speed processing there is no easy method to gain 
access to the actual intermediate steps that are taking 
place in the long chain of signal processing functions. 

A solution to this problem of student access to the 
intermediate steps is to simulate the signal processing on a 
personal computer. This simulation can be designed so as to 
show any step in the process that is desired and to allow 


the student to vary the inputs at will to compare different 


algorithms. The simulation also calculates the performance 
characteristics particular to the process under 
consideration to allow comparisons between different 
algorithms. It is often easier for a student to comprehen 
visual outputs rather than analyze a table of performance 
characteristics, so the simulation will display the RSP 


results graphically. 


B. FUNDAMENTALS OF RADAR SIGNAL PROCESSING 
Older non coherent radars could detect only targets 

whose input signal to the radar was larger than the 
Surrounding background interference. Smaller target signals 
hidden in the ground return, rain clutter and oem 
interference were lost to the system since it was not 
possible to isolate the target signal from these 
interferences. This concept is shown in Figure 1 where only 
Target 1 could be detected by a non coherent radar. 
However, modern radars utilize the doppler effect to isolate 
the weak Signal of a moving target from much larger 
background interference. This allows the detection of 
signals in interference that the older radars could not 
accomplish. Thus Targets 2 and 3 from Figure 1 could be 
detected by a modern radar. Implementation of this concept 
is shown in Figure 2 and its sub components are as follows 

® The programmable signal processor (PSP) is the first 


computer used in the signal processing chain. It is 
designed as a high speed, computationally intensive tool 


Pomprecess the Signal Erom the receiver input (Rx). 


@ The A/D converter samples the signal at a rate on the 
order of 1 MHz. 


® The adaptive nulling is used to cancel any interference 
coming in through the sidelobes of the radar’s antenna. 


@® Pulse compression is used to compress the long 
transmitted pulses in order to have finer range 
resolution. This stage develops a processing gain equal 
to the pulse compression ratio (PCR) and thus builds up 
the target signal against background interference 


® The clutter canceler significantly reduces ground 
Sluttem addmer len stamrondt, TOmsaeor Clutter using the 
moving target indicator (MTI) concept 


® The doppler processor takes the spectrum of the incoming 
Signal by performing fast Fourier transforms (FFT) on 
the signal samples. This process isolates the small 
target Signal, which has undergone a doppler shift due 
to target velocity, from any Stationary interference. 


@® The constant false alarm rate (CFAR) detector examines 
increases in signal amplitude above the average 
background interference in order to detect targets. 


® The signal processing up to this point has been very 
computationally intensive and must be done ina high 
speed, optimized processor known as a programmable 
Signal processor (PSP) in order to detect targets in 
real time. 


® Data rates after the CFAR detector are much reduced and 

the processing is then done in a slower computer known 
as radar data processor (RDP). This processing is 
performed across several coherent pulse intervals 
(CPI’'s), and the end result is the target report which 
normally consists of the target’s range, azimuth and 
meloOcity. 
These signal processing steps are "invisible" to the 

student as there is no practical way to access these stages 

in either of the computers. This is the area which most 


feeder simulation in order for students to become familiar 


Wael it. 


Target 1 
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Figure 1 : Targets in clutter for a non MTI analog radar 
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Figure 2 : Modern radar block diagram 


C. TRADITIONAL METHODS AND SIMULATION 

The overall subject of RSP consists of the underlying 
theory and implementation of digital signal processing 
algorithms. The particular problems inherent in the 
teaching of RSP in the classroom that can be overcome by the 
use of Simulation are explained in this section. 

1. Textbooks 

The textbooks available in this area do a good job 
in explaining the basic principles of the subject. However, 
by necessity textbooks are limited to discussing only a few 
cases on the topic at hand due to limited text space. Even 
these few cases will require extensive computations which 
can not be shown in their entirety. 

Simulation relieves the student of this 
computational burden. This leaves the student free to 
explore any aspects of the process which interests him. The 
student can experiment by changing the inputs or parameters 
of the problem and see their effects visually. Practical 
demonstrations on the application of the theory aids and 
motivates the student in comprehending the underlying 
mermciples of the process. 

For example, in the section on pulse compression the 
Student can easily find that not all signals will compress 
in a useable fashion, especially under doppler conditions. 


All he has to do is choose a signal and the doppler shift it 


is to undergo, and the simulation will give the vous outa 
that particular case. The result is shown in graphical form 
and this brings the point home immediately. 

2. Laboratory Hardware 

In actual radar systems, as well as those set up in 
laboratories, the signal processing functions are performed 
by dedicated computers such as the PSP and the RDP. These 
computers pass the data to each other and sometimes send 
Signals to other components of the radar. However, there is 
no easy way to access the outputs of each separate signal 
processing function. The intermediate results are almost 
completely inaccessible, even in a laboratory environment. 
Thus the most important parts of the RSP processes are 
hidden from the student. 

This problem is also solved by using simulation 
programs to simulate each RSP process separately. This 
allows the student to separately examine each process, vary 
its inputs and examine the changes which occur within that 
process. He can then combine them for end to end 
Simulation. 

3. Numerical Computations 

It is a fact of life that normal teaching methods 
lack the ability to study realistic radar problems. The 
following illustration on the detection of moving targets at 


a Single azimuth angle explains why. A typical computation 


would consist of perhaps 5 coherent processing intervals 
(CPI) for a MTI radar. Each CPI would typically consist of 
1000 range bins for each of the 16 pulses constituting a 
CPI. Thus data for one azimuth angle consists of 5 x 1000 x 
16 = 80,000 complex numbers, and each component of a complex 
number is typically represented by 12 bits. The doppler 
processing alone would require the computing of 5,000 
Sixteen point FFTs resulting in 80,000 range doppler cells. 
This doppler processing would be followed by a chain of 
algorithms and thus the processing for just one azimuth 
angle is well beyond a student’s capability. 

However, this signal processing can be readily 
Simulated on a personal computer Since Simulation does not 
need to be done in real time aS in an actual radar. [It 
takes a few more seconds but produces results of the same 
accuracy as the real world system. 

4. Simulation of RSP 

There 1s one key point about the Simulation of RSP 
algorithms which must be noted. The simulation will be just 
as accurate as the real system would be given the same 
inputs. The only difference between the real system and 
this simulation is that the generation of the input signals 
is done by the computer. All other algorithms and processes 
are performed in the same manner as they would be done in 


the PSP or RSP, albeit that the personal computer used for 


these simulation programs is many orders of magnitude slower 


than the PSP. 


D. REVIEW OF LITERATURE 
Radar signal processing is a specialized application 
within the signal processing area. While RSP is an integral 
and critical part of a radar system, there are very few 
texts or works devoted solely to RSP due in part to the 
difficulty of discussing it uSing the traditional teaciam 
methods. In general the best sources are either texts on 
general radar applications or those on signal processing. 
In terms of the topics discussed in this thesis, the 
best source works were found to be the following: 
® DFT as a filter bank [Ref. 1: Part V], [Ref. 3: chap 
14], (Ref. 4: chap 7], {Ref. 6:p 119 to 125], {Rema 
chap 7], (Reb. 36: bp 1233. 6o S27 
® Weighting schemes for a DFT filter bank [Ref. 3: chap 
10.8], (Ref. 6:p 426 to 427], {Ref. 7: chap 7 ange 
[Ref. 9: Appendix D]. 


® Pulse Compression [Ref. 1:Part IV], [Ref. 2:chap 2], 
(Ref. 3:chap 4 and 10], [Ref. 4:chap 8], [Ref. 5:chap7]. 


The above list is by no means complete, and better 
sources for individual topics may exist. The thesis in 
itself should provide a good introduction to these topics 
for a student in RSP. In particular the simulation pregmaie 
will allow him to manipulate the RSP algorithms to gain a 


fuller understanding of the processes involved. 


tLe: 


el. THE DISCRETE FOURIER TRANSFORM AS A FILTER BANK 


Modern coherent radars discriminate between targets and 
clutter on the basis of target motion. Moving targets 
receive a doppler frequency shift that causes the frequency 
of the received target return to be different from the 
Stationary clutter return. Signal processing methods that 
take advantage of the frequency shift for moving target 
indicator (MTI) radars are the delay line canceler and the 
discrete Fourier transform (DFT) filter bank. 

The delay line canceler for MTI radar 1s reviewed and 
the advantages of the DFT filter over it are highlighted. 
The development of the DFT as a filter bank for radar 


applications is addressed more thoroughly. 


A. DELAY LINE CANCELER VERSUS DFT FILTER BANK 

The delay line canceler used in MTI radars is a time 
domain filter. It rejects the clutter at low frequencies 
and passes the target signal at higher frequencies. It is 
essentially a band reject filter with the notch placed at 
dominant clutter frequencies. [Ref. 6:p107] One simple way 
of doing this is by storing the input value for a period of 
time equal to the PRI and then subtracting the stored value 
from the new value. This removes the stationary clutter at 


low doppler frequencies. However, it also cancels any 


ibal 


Signals at frequencies that are close to multiples of the 
PRF. Figure 3 shows the spectrum of a radar return with 
ground clutter, rain clutter and targets ae Gnesranee 
interval. The radar can not separate out the return Bbyewa 
of the doppler frequencies and thus processes all of the 
frequencies present in this interval from zero to the PRF. 
The ground and rain clutter are present at all times for 
this range and thus the constant false alarm rate (CFAR) 
mechanism raises the detection threshold to a level that 
requires a greater target magnitude response than the 
average clutter level. The detection threshold for this 
Situation is shown in Figure 3. Only Target 1 could be 
detected as it is the only target that exceeds the 
background clutter level. 

Figure 4 shows the spectrum at the output of a delay 
line canceler where the ground clutter has been attenuated. 
The radar still processes all of the frequencies present in 
this interval from zero to the PRF, but the magnitude of the 
stationary clutter has been attenuated so the CFAR detection 
level is much lower than in Figure 3. However, the delay 
line canceler has no effect on the rain clutter, and Target 
2 is still masked by the rain clutter since the CFAR 
detection level is now dictated by the rain clutter 
magnitude. However, if a DFT filter bank were used, then 
the rain clutter can be isolated from target 2. This is 


shown in the spectrum of the return for an eight filter DFT 


a 


bank at Figure 5. Now the rain clutter is ina separate 


meeter Erom the return of Target 2. Thus Target 2 could now 
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Figure 3 : True spectrum of a radar return showing the 
detection threshold set by CFAR 


be detected as the CFAR is set individually for each filter 
and Target 2 is significantly greater than the thermal noise 
in its filter. The main advantages of a DFT filter bank 
over a delay line canceler are: [Ref. 6:p107 to 117], 

Mer 5:p7 to 33) and (Ref .12] 

@® The spectrum of the return is divided into several 
different filters to better separate clutter from 
targets. 

@® Also, it is possible to change the filter width by 
varying the PRF or N so that any targets that appear 
near clutter frequencies for one PRF can be moved to a 
filter that has very little clutter in it for a second 


PRF. This will be demonstrated more clearly in the 
Simulation program in Chapter III. 
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Figure 4: Spectrum of a radar return at the output of a delay 
line canceler. The ground return clutter and any returns at 
frequencies that are near multiples of the PRF have been 
canceled. » 
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Figure 5 : Spectrum of radar return after a DFT filter bank. 
The ground clutter is in filters 0 and 1, while the rain clutte 
is in filters 3 and 4. Target 1 is in filter 2 and Target 2 is 
in fa Leenes. 
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@® Receiver noise is limited to one filter width and thus 
the signal to noise (SNR) level is raised for a target 
meOEUIN . 
@® The possibility exists to measure the target's doppler 
frequency and thus compute its radial velocity. 
B. DFT FILTER BANK FOR COMPLEX SINUSOIDAL INPUT 

The development of the DFT as a filter bank in radar 
applications is presented for a complex sinusoidal input and 
for a real Sinusoidal input. A complex sinusoidal input 
would result from the use of quadrature sampling of an input 
signal whereas a real input would be from a single channel 
sampling process. Figures 6 and 7 show how these inputs are 
formed. The detector in these figures consists of a mixer 
followed by a low pass filter. The main advantages of 
quadrature sampling over single channel sampling are: 


@® Positive and negative doppler frequencies can be 
distinguished. 


® The blind phase problem can be avoided in quadrature 
sampling. Blind phases may occur in single channel 
sampling when the sampling interval is such that it 
coincides with the zero crossing of the incoming signal. 
For example, when ¢=90°, cos(¢?) is a very small number 
for single channel sampling. In quadrature sampling 
these blind phases are eliminated by the presence of 
sin(¢@), which will be high when cos(¢@) is low. [Ref. 
Peel S.7)], [Ref. 5:p7 to 40). 


The development of the DFT for a quadrature sampled signal 
PorOws . 
1. The Input Signal 
If the input signal is a complex sinusoid it can be 


represented by 


Ls 






Input Signal 


A cos(ot+o) 7 


Figure 6 : Quadrature Sampling which results in a complex 
input signal 







Input Signal 


A cos(wt+ 6) 


Figure 7 : Single Channel Sampling which results ina real 
input sinusoid 
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S102) = OOS 2a) = 4) Soo eyegie  — eau (1) 


Once it has undergone quadrature sampling it can be 
represented as a discrete signal by noting that x = fNT and 


m@e= nt so that 


® x is the discrete frequency. 
® N is the total number of samples. 
@®n is the sub sample number, n = 0O:N-1, nel. 


@® PRF is the pulse repetition frequency. This is the 
sampling frequency for radar applications. 


@® nT is the discrete sample time. 
2. Discrete Fourier Transform Development 
The DFT for a finite length sequence s(n) is as 


Eollows 


N-1 ~7 2nnk 
F(k) = }> s(n) e” *% (3) 
n=0 


e@eseituting for s(n) from Equation 2 leads to the DFT 


response for a complex sinusoidal input 


N-1 j 2nn(x-k) N-1 
F(k) =-Ye N = G2, (4) 


n=0 n=0 


ay 


(SF 


where 8 = AE 


Then using the mathematical relation for the summationmecguee 
Equation 5 can be rewritten as 


finite geometric series, 





1-e- 
eae 
= © nyo -e 
,) 8 8 
-j= = a3) 6 
e * e *-e + (6) 
NO 
jw sin(=>) 
= -3 2 * ) 
Sint) 
ze 


Magnitude Response of the DFT to a Complex Input 


ce 
In the detection of targets only the magnitude at 
and thus phase can be 


the £1liter OurpUE 1s Considered, 
The magnitude of the output of the DFT will be 


ignored. 
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The frequency response of the DFT to a quadrature sampled 
input Sinusoid is shown in Figure 8. This plot is of the 
magnitude response of Equation 8 for k = 0 over a frequency 
range of several PRFs. The magnitude response pattern of 
F(k) is repetitive in nature and repeats itself every PRF 
because the PRF is the sampling frequency. This then leads 
to the conclusion that F(k)=F(k+N). This is proven as 


follows: 


IF (k)| = |FCksN) | ) 


N 
= | Sees hess) ue 
Shon ( DS en) 
N 
a S akaat is (ose 
Pt ( BAX=K) (12) 
; N 
= EWS) OED (giles) 


The same response for k = 0 is shown in Figure 9 
over one PRF from the SINGPLOT option of the DFTBANK 
Simulation program found in Chapter III. This pattern is 
Gonsidered to be the 0" filter (ie. filter # k=0) of a DFT 
filter bank. If the magnitude response of this filter 
crosses the detection threshold, then a target is reported 


as being within the frequency range of the mainlobe of the 


io 


magnitude response for that filter.» Showld a target iim 
Sidelobes of the filter exceed the threshold value, then the 
radar will assume it is coming from within the mainlobe and 
thus a false target frequency will be declared. 

The null to null frequency width of the filter isis 
on Figure 9, as is the 3 dB frequency width. The equations 


to determine these widths are 


AN oe 


nuia 


= Null to Null MAINLOBE WIDTH = 





Mgega (14) 
N 


and 


Af,qp = 3 dB MAIN LOBE WIDTH = 0.89=—=. (15) 


By varying k from 0 to N-1, kelI the DFT filtereeage 
is formed from the mainlobes of the N filters. This is 
Shown in Figure 10 from, the MULTMAIN option of the DFTBANK 
Simulation program. The straddle loss is defined as the 
loss in magnitude response between the peak response of the 
filter and the response of the mainlobe when it crosses the 
mainlobe of the adjacent mainlobe. The frequency range for 
mainlobe detection of a filter is defined as being from one 


adjacent crossover point to another. This is shown in 


20 
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Figure 8 : F(k) Magnitude Response for k = 0 over several 
PRFs. The pattern repeats itself every PRF. 
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Figure 10 : F(k) Magnitude Response for k = 0 uSing the 
SINGPLOT option, PRF = 2000 Hz and N =8. 
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Figure 9 : DFT filter bank showing filter mainlobe 


crossovers for PRF = 2000 Hz and N=8 from MULTMAIN option 


DFTBANK simulation. 


Zo 
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Bagure 10. 


Thus the crossover point defines the worst 


Af = 


response a target would have while within a filter bank. 
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DFT FILTER BANKS FOR REAL SINUSOIDAL INPUT 
ae? 


The Input Signal 


If the input Signal iS a real sinusoid then it can 
be represented by 


Coe ie 


cos (2nft) 


a8 [esente i @ ~Jente) 
2 


(ie 
The sampled signal can be represented in discrete form by 


(18) 
2. Discrete Fourier Transform Development 


The form of the DFT for the real finite length 
sequence s(n) 


is the same as in the complex case. 
Substituting Equation 18 into Equation 3 leads to 


Thus 
N-1 l < 2nn(x-k) a 2nmn(x+k) 

F(k) = Dasube oY +e Nie Des Be (19) 
n=0 


Ze 


Geyer tne Guality inherent in Equation 19 it is evident that 
F(k)=F(-k). Using the mathematical relation for the 
convergence of a finite geometric series in the same fashion 


as for the complex sinusoid results in 


ois .  » NO 
a = neay sin(—) <5 (RE) sin(—) 
F(k) ==le a ee he ¢*——_—_], (20) 
¢ sin(—) sin(®) 
2 2 
where @ = 2% — (21) 
and 6 = oa BE 2) 


3. Magnitude of the DFT Response to a Real Input 


Rewriting Equation 20 leads to: 


=~7(N-1) 2 ee 
Bik) = ee - 5 
° sin(S) 
(28)) 
Ss Uene yee sin(%®) 
1p lS 2 
ee) 
2 


Therefore the magnitude response is 
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s eh 2 
PES sce eal aaa: 
Sint) 
(24) 
5 tw-1) £o-8) sin( 32) 
J |. 
ce) 
sin(—) 


term can be reduced by 


The remaining exponential 
@ from Equations 21 and 22 so 


substituting in for @ and 


x=) 28 ¢ (04%) ~(x-k)) 


~>(N-1) 
—— ee oe 
N 
= @ 
onk (25) 
= @J2tKy N 
x 2K 
= Be N 


Substituting Equations 25 into Equation 24 thus gives a 


magnitude response of 


i ae sabes) | 
sin (ZAXt*) | 
N 
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The frequency response of the DFT for a single channel 
sampled input sinusoid is shown in Figure 11. This plot is 
of the magnitude response of Equation 26 for k = 1 using the 
SINGPLOT option from the DFTBANK simulation program at 
Chapter III. For k=1 there are two equal magnitude 
mainlobes spaced evenly about a centre frequency that is a 
multiple of the PRF. The obvious difference between the 
complex and real input data cases is that the former has a 
Single main lobe and the real sinusoid has two main lobes 
that are equally offset from the PRF. In Figure 12 for k=-1 
the magnitude response of the DFT is identical to that of 
Figure 11. This is as expected given the duality of 
@eam@atron 19 such that F(k)=F(-k). Therefore it is not 
possible to distinguish positive and negative frequencies of 
the same magnitude since they will have an identical 


magnitude response from the DFT filter bank. 
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Figure 11 : F(k) magnitude response of a real input sinusoid f 
k=1 using the SINGPLOT option, PRF = 2000 Hz and N=8. 
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Figure 12 : F(k) magnitude response of a real input sinusoid 
for k=-1 using the SINGPLOT option, PRF = 2000 Hz and N=ee 
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Hite oe SIMULATION OF A DFT FILTER BANK 


This chapter describes the methodology of the simulation 
of a DFT filter bank response to a real or complex sinusoid, 
presents an overview of the program structure, describes how 
to use the program and shows how to interpret the results. 

The DFT filter bank simulation program implements a 
filter bank showing the magnitude response of a real or 
complex sinusoidal input for all unambiguous frequencies. 
The following parameters can be varied in the simulation; 


® N = the number of points for which the DFT is to be 
calculated, 


@ PRF = the pulse repetition frequency of the radar, 
@®k = the filter designator (if required), and 


meee — Ene doppler frequency of the target (if 
required). 


By controlling these parameters the student can 
Semeeruct DFT filter banks for inputs that are real or 
complex sinusoids. Some of the situations that can be 
addressed are: 

@® The response of a filter bank to a target appearing at a 
feven doppler frequency (£,,) can be assessed. If the 
target’s doppler frequency 1S an ambiguous frequency of 
the filter bank, then the alias target frequency (f,) 


report is given. The alias target frequency is 
seveloped from 


ae 


f, = £.,, modula Gla (277) 


@® If the target and clutter are in the same filter ae 
given PRF and N, then the filter bank can be changeduby 
varying either N or the PRF to move the target and 
Clutter into separate filters. 

® The response of a DFT filter bank using Single Ghanem 
or quadrature sampling to positive or negative doppler 
targets can be observed. 

® There are numerous other possibilities left to the 
user’sS initiative. 

A. PURPOSE OF THE SIMULATION PROGRAM 
The equations implemented by the subprograms are fora 
real input Sinusoid from single channel sampling anda 


complex sinusoid from quadrature sampling as seen in Figures 


6 and 7. The frequency range of the filter bank is for 


(28) 


\f| < Prom 
2 


which covers all of the unambiguous frequencies for that 
filter bank. The frequency is stepped in increments of 
0.005*PRF which gives a smooth plot while maintaining a 


reasonable processing speed on a DOS 386 machine. 


B. PROGRAM STRUCTURE 

The overall program structure can be seen in Figure 13. 
The structure is that of a main program (DFTBANK) which 
calls various subprograms (SINGLOBE, MULTMAIN, ... etc) as 


required. There is an error message (ERRORMSG) which 


Zo 


detects an improper data entry, notifies the operator and 
restarts the process. It is possible to exit the simulation 
from either the main or subprograms. This structure was 
used since it allows maximum operator control and input. It 


is a robust format that 1s data entry error tolerant. The 
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Figure 13 : Overall program structure for DFT bank simulation 


other advantage of this format is that if an operator 
desires to change one of the subprograms or add another 
subprogram the system can accommodate this. Any errors made 
in this process will affect only the subprogram in question. 
The rest of the simulation will remain operational while the 


new addition is being debugged. 
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OPERATING INSTRUCTIONS AND GUIDE 
1. Initial Setup 


To use this simulation the following MATLAB.m files 


must be copied into the user sub directory. Copies of these 


programs are available from Dr GS Gill at USNPGS. 


DFTBANK.m is the simulation control program. 


SINGLOBE.m 1s the simulation subprogram responsible for 
displaying the magnitude response of a single filter 
across all unambiguous frequencies in the filter bank. 


MULTMAIN.m is the simulation subprogram responsible for 
displaying the mainlobe magnitude response of all 
filters in the filter bank for all unambiguous 
frequencies. 


MULTLOBE.m is the simulation subprogram responsible for 
displaying the full magnitude response for all filters 
in the filter bank for all unambiguous frequencies. 


TGTPLOT.m is the simulation subprogram responsible for 
calculating and displaying the filter response for a 
target at a given doppler frequency. It can accept 
ambiguous target frequencies and the filter response 
will be for the alias frequency. The filter that this 
alias frequency lies in can result in a positive doppler 
being identified as negative or vice versa. 


ERRORMSG.m informs the operator when an entry error has 
been made and restarts DFTBANK after clearing the system 
of all variables. 


This simulation works on a DOS 286 machine, but a 386 68 


higher is preferable for increased speed. To start the 


Simulation type DFTBANK [ENTER] at the MATLAB prompt. 


2. Printing Graphical Outputs 


The assumption has been made that the PRINT SCREEN 


key will be used to print out plots. To use this key the 
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DOS graphics command must have been issued at the DOS prompt 
before any printing is attempted; C:\>> graphics [ENTER] is 
the necessary command. 

Should it be preferable to create meta files for each 
plot, the operator has two options: 


® Modify all plot options in the subprograms by adding a 
meta command and file name after the plot command, or 


® Press the Ctrl-C keys to exit the simulation once the 
plot is on the screen. This will return him to the 
MATLAB prompt and he can issue the meta command from 
there. To restart the Simulation type DFTBANK [ENTER]. 
3. Program Options 
After the DFTBANK command has been issued to start 
the simulation the first screen seen by the operator is at 
Figure 14. The second screen the operator will see is at 
Figure 15. At this point the operator must choose whether 
he wishes to process a real or complex sinusoid for the 
Simulation. Then he has to choose the type of system 
response that he wants to manipulate. Once this has been 
done the main program transfers control to the subprogram 
selected in Figure 16., All of the screens for the 
subprograms are very Similar, the only difference being 
eer "K"sor ME," are required inputs. 
The choice for a dB or ratio plot determines whether 
the magnitude scale on the y-axis will be plotted for ratio 


values or in a decibel format. The dB format is preferable 


when determining the sidelobe level response. For a more 
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This program will allow you to create Doppler Filter 
Banks for either real or complex input sinusoids. Real 
Sinusoids are the result of single channel sampling while 
complex sinusoids are from quadrature sampling. Single 
channel sampling means that it will be impossible to 
differentiate between positive or negative doppler. 
Quadrature sampling allows positive or negative doppler 
to be distinguished. 





The options are to see the full response of a single 
filter over the full spectrum of unambiguous frequencies, 
to see the mainlobe response of a DFT bank of filters, to 
see the full response of a DFT bank of filters (including 
| sidelobes) or to see the magnitude response of each 
filter in a DFT bank to a target occurring at a 
particular frequency. 





TOUCH ANY KEY TO sGONT INGE 


Figure 14: First Screen of DFTBANK main control program 


You must choose whether to process a real sinusoid 
from single channel sampling or a complex sinusoid from 
quadrature sampling. 


FOR A REAL SINUSOID, choose 1 
FOR A COMPLEX SINUSOID, choose 2 
ENTER YOUR CHOICE NOW, SINUSOID = 


You must choose which type of response you wish to see 
To see a Single filter response, ENTER 1 

To see a filter bank with main lobe response,ENTER 2 

To see a filter bank with side lobe response, ENTER 3 

To see target values calculated for each filter,ENTER 4 
TO EXIT THIS PROGRAM, ENTER 0 

ENTER YOUR CHOICE NOW 

CHOLGE= 





Figure 15: Second Screen of DFTBANK main control program 
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general overview, the ratio format is better. The actual 
inputs for the variables described in Figure 16 are given on 


the next screen at Figure 17. 


YOU MUST NOW ENTER SOME OF THE NECESSARY PARAMETERS. 

These parameters are the number ’N’ of DFT points the 
filter is to calculate. This is also the number of 
filters there will be in the filter bank. 

Next you must choose the PRF you wish to consider. 

Then you must choose the ’k’ offset for the DFT to 
determine which filter in the DFT filter bank you will 
select. The choice of ’k’ is an integer value usually 
from 0 to N-1, but you are not limited to this selection 
if you desire to see the periodic nature of the DFT. 

Finally you must choose whether you wish the display 
of the magnitude response to the DFT to be in normalized 
Semor ratio form. 


TOUCH ANY KEY TO CONTINUE 





Figure 16: First Screen of a subprogram 


Number of DFT points you wish is, 
The PRF in Hz you wish is, PRF = 


The Filter K # you wish to use is, 
mer a G8 plot choose 1, for a ratio plot choose 2, 





Figure 17: Second Screen of a subprogram input request 

Once all of the graphical outputs and results have 
been shown, each subprogram will terminate with the screen 
at Figure 18 where the operator can choose to do the same 
Subprogram for the same type of input to the DFT (the 
SINUSOID variable remains constant), return to the main 
program (DFTBANK) to try a different simulation or exit the 


Simulation. 


ene 


Any data entry error that the inherent MATLAB error 
routines do not recognize will be detected by the simulation 
and the error message at Figure 19 will be displayed. All 
variables will be cleared and the main program DFTBANK will 


be restarted. 


—_ — ES — = 


YOU MUST NOW CHOOSE TO EITHER DO THIS CHOUCE AGA as 
RETURN TO THE MAIN MENU OR EXIT THE SIMULATION 






DO THE SINGLE LOBE AGAIN, CHOOGSE ai 
TO RETURN TO THE MAIN MENU, CHOOS EaZ 
TO BALE Tob Siva eerie CHOOSE 0 






Figure 18 : Final Screen of a subprogram 


YOUR CHOICE DOES NOT CORRESPOND TO THOSE AVAILABLE. 


THE PROGRAM WILL BE RESTARTED FOR YOU. 
REMEMBER YOU MUST CHOOSE O TO EXIT. 


PRESS ENIERK TO CONTINUE]. 





Figure 19 : Error message subprogram screen 


D. INTERPRETATION OF RESULTS 

The output of the simulation is mostly graphical in 
nature Since graphical displays are readily assimilated and 
provide a more intuitive understanding of the process. 
Results from each of the subprograms are presented. These 
plots are analyzed to demonstrate some of the capabilities 
of the simulation program and to provide worked examples as 
a reference guide. In each case the data inputs used are 


provided and the selection process described. 
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1. Single Filter Response 

The SINGLOBE subprogram displays the magnitude 
response of a single DFT filter across all unambiguous 
Prequencies in the filter bank. These unambiguous 
frequencies are for |f |s PRF/2. While it is the mainlobe 
Magnitude response of the filter that defines the frequency 
coverage of the filter, there will be some response to any 
frequency present in the filter bank. If the response for a 
frequency were from outside of the mainlobe magnitude 
response range and it exceeded the detection threshold set 
for that filter, an additional target report would be 
generated. 

The operator can vary the PRF and N in the 
Simulation to achieve a filter response with an acceptable 
filter width, frequency resolution and doppler frequency 
coverage. In the following example the inputs to the 
Simulation were: 


SeotNUSOID = 1 for a real input Sinusoid to be processed, 


CHOICE = 1 for a single filter response, 


@ PREF = 2000 Hz, 
@® k = 1 for the F(1) magnitude response, and 
@® dB ratio = 2 for a ratio plot. 
All of these inputs will stay the same for all 


subsequent simulations in this section unless noted 


BS 


otherwise. The plot resulting from these simulation choices 
is at Figure 20. The twin mainlobe response for a real 


input sinusoid is clearly evident. Had dB ratio = 1 been 


Single Lebe of oe B Point OFT Filter Bonk, PRE =» 2000 He 








FCk) response for «amt for fo red! (nput ainusoid 


Wagnilude 


Figure 20 : Magnitude response of a real input sinusoid for 
F(k) in a DFT filter bank where N=8, PRF = 2 KHz and k=1 


selected than the dB ratio magnitude response at Figure 21 
would have been plotted. Selecting k=-1 would result in 
Figure 22. This shows why positive and negative doppler can 
not be distinguished by a real input DFT filter bank as 
F(k)=F(-k) so that targets at the same positive or negative 
doppler frequency will have the same response. 

Selecting k=7 would result in Figure 23, which shows 
the recurring nature of the DFT as F(k)=F(k+N) with respect 
to Figure 22 where kK=-1. Had SINUSOID=2 for a complex 
Sinusoid from quadrature sampling been chosen then the graph 
at Figure 24 would have been plotted. The differences 


between real and complex input sinusoids are evident. 
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Figure 21 : Magnitude response in dB of a real input sinusoid 
Bort F(k) ina DFT filter bank where N=8, PRF = 2 KHz and k=1 
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Figure 22 : Magnitude response of a real input sinusoid for 
F(k) in a DFT filter bank where N=8, PRF = 2 KHz and k=-1. 
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2. Mainlobe Filter Response 

The MULTMAIN subprogram displays the mainlobe 
magnitude response of all filters in the filter bank forme” 
unambiguous frequencies. The choices are the same as for 
SINGLOBE, with the exception of "k" since the mainlobe 
magnitude responses for k=0 to N-1, k ¢€ I are all plotted 
automatically. 

In the following example the inputs to the 
Simulation were 

® SINUSOID = 2 for a complex input sinusoid, 


@® CHOICE = 2 for a filter bank with mainlobe response, 


® PRF = 2000 Hz, and 

® dB raeio = 25for a rare ler. 
All of these inputs will stay the same for all subsequent 
Simulations unless noted otherwise. The plot resulting from 
these simulation choices is at Figure 25. The crossover 
filter width, straddle loss and other performance 
measurements can be noted from it. Should SINUSOID = i have 
been chosen then the plot at Figure 26 would have resulted. 
The different filter bank peak responses to a real input 
Sinusoid for the various k inputs of Equation 34 can be 
seen. Therefore the detection thresholds for the different 


filters must vary to account fon Enis yore 
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Figure 23 : Magnitude response of a real input sinusoid for 
F(k) in a DFT filter bank where N=8, PRF=2 KHz and k=7 
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Figure 24 : Magnitude response of a complex input sinusoid for 
F(k) in a DFT filter bank where N=8, PRF=2 KHz and k=1 
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3. DFT Filter Response With Sidelobes 
The MULTLOBE subprogram is the same as the MULTMAIN 
one, except that the full response for each filter is 
calculated and plotted, not just the mainlobes. 
In the following example the inputs to the 
Simulation were 
@® SINUSOID = 2 for a complex input sinusoid, 
@® CHOICE = 3 for a filter bank with side lobe response, 
@N = 8, 
@® PREF = 2000 Hz, and 
® dB ratio = 2 for a ratio plot. 
All of these inputs will stay the same for all subsequent 
Simulations unless noted otherwise. The plot resulting from 
these simulation choices is at Figure 27. This demonstrates 
where the various filters have a mainlobe crossover and the 
peak sidelobe levels. From it the choice of where to set 
the detection threshold can be based on the following; 
® to avoid declaring targets in more than one mainlobe of 
a filter, the threshold must be set higher than the 
MNdanlobe Crossover poi: aor 
® to avoid declaring targets from the sidelobe of a 
filter, the threshold must be set higher than the 
highest sidelobe. 
This gives the operator a good feel for the dynamic range 


over which he can detect targets that vary in signal 


Strength and not have excessive detection threshold 
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Figure 25 : Magnitude response of a complex input sinusoid for 
F(k) in a DFT filter bank where N=8, PRF=2 KHz and k=0 to N-l, 
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Figure 26 : Magnitude response of a real input sinusoid for 
F(k) in a DFT filter bank where N=8, PRF=2 KHz and k=0 to N-1, 
kel 
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crossings and false doppler reports. A measurement of this 
nature 1S easier to do if the dB ratio = 2 for a Gags plomeueeee 
been selected as in Figure 28. 

4. Filter Bank Response to a Specific Target Frequency 

The TGTPLOT subprogram is the same as the MULTLOBE 

Subpogram except that 1t will require an additional input 
for the doppler frequency the target or clutter 1S COWamee ee 
at and that all plots are automatically of a dB nature. If 
the doppler frequency given does not satisfy the Nyquist 
criteria then an alias frequency (f,) will be calculated 


from 


f,=f 


: op MOCULOW ae (29) 
This alias frequency is then processed by the filter bank. 
The filter that the alias frequency lies in can result ina 
positive doppler being identified as negative or vice versa. 
In the following example the inputs to the 

Simulation were: 

@ SINUSOID = 2 for a complex input sinusoid, 

® CHOICE = 4 £oOr a targee plore 

@® N = 8 

® PREF = 2000 Hz. “and 


@ > £2— 4725.82 


gt 
The plot resulting from these simulation choices is at 


Figure 29. The plot of Figure 30 shows how a [= 95250 


42 


that does not satisfy the Nyquist criteria can result ina 


f, that gives misleading results. 
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Figure 27 : Magnitude response of a complex input sinusoid 
for F(k) ina DFT filter bank where N=8, PRF=2 KHz and k=0 to 
N-1, keI showing all sidelobes 
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Figure 28 : Magnitude response in dB of a complex input 


Sinusoid for F(k) ina DFT filter bank where N=8, PRF=2 KHz 
and k=0 to N-1, keI showing all sidelobes 
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Figure 29 : Target response for each filter of a complex 
Pijede sinusoid DFT filter bank where N=8, PRF=2 KHz and f,,,=425 
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Figure 30 : Target response for each filter of a complex 
Miput Sinusoid DFT filter bank where N=8, PRF=2 KHz and 
meee 50 Hz. Note that £,=-750 Hz so anegative doppler 
target 1S incorrectly declared. 
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IV. WEIGHTED DFT FILTER BANKS 

As seen in Chapters II and III the sidelobe response of 
a DFT filter is significant enough to cause problems in the 
detection of weak targets while maintaining a low false 
alarm rate. This problem results from the fact that a 
Strong return from a single target could have enough of a 
magnitude response in the sidelobes of the other filters to 
exceed the detection threshold for those filters. Thus 
false target reports are generated, and the radar sees 
multiple targets when only one exists. If this problem is 
solved by raising the detection threshold to the point where 
the sidelobe response of a large target no longer causes 
false target declarations, then the radar has been 
desensitized to targets with smaller magnitude responses. 

What 1S required is a method to lower the sidelobe 
response level of a DFT filter to an acceptable level. One 
way of doing this is through the use of windows (or 
weighting). This technique is in widespread use for 
communications systems and for other parts of radar signal 
processing such as reducing the time sidelobes in pulse 
compression. This chapter concentrates on some of the 
Simpler DFT filter weighting schemes and describes the 
window functions that produce them. This chapter 1s Tinimiees 


to non adaptive weighting schemes in common use, such as: 
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and 


Rectangular window, where the inputs are uniformly 
weighted. The frequency domain transfer of the time 
domain weighting sequence is known as the Dirichlet 
meommecl, PD. (K) . 


Bartlett (or Triangular) window where the inputs are 
weighted by a triangular function. 


Hanning (Von Hann, or raised cosine) window where the 
inputs are weighted by a raised cosine function such 
that the magnitude of the response at either end of the 
window is at zero. 


Hamming window which manipulates the Dirichlet kernel to 
further reduce the sidelobes immediately adjacent to the 
main lobe. 

Blackman window is an extension of the Hamming window in 
that 1t manipulates more terms of the Dirichlet kernel 
for greater Sidelobe reduction. [Ref. 9:p566 to 574], 
feet. 14:p98 to 105] 

WINDOWING 


The DFT accepts an input sequence of N samples, s(n), 


processes them according to 


N-1 _7 2mnk 
F(k) =} s(n) e™ * (30) 
n=0 


A windowing function, w(n), weights each input sample by a 


specific amount. The amount each input sample is weighted 


is a function of the windowing scheme used. The new input 


sequence to the DFT is s,(n), and it is formed by 


The 


s(n) =w(n)s(n) for n=1-N, nel. (31) 


W 


main benefit of the windowing scheme is a reduced 


sidelobe magnitude response; the chief disadvantage of a 


windowing scheme is that the mainlobe of the filter is 


4‘7 


broadened. The penalty of mainlobe broadening due to 
windowing is that the filter width is increased; therefore, 
the ability to resolve targets and clutter close together in 
frequency 1s degraded. Also, there is a greater incidence 
of a single target being reported in adjacent filters if N 
1s held constant. 

The method to decrease the filter width to an acceptable 
level is by increasing N or by decreasing the PRF. Thus the 
search for better windowing functions has concentrated on 
finding windows that greatly reduce the sidelobe magnitude 
response, have as narrow a mainlobe as possible and are easy 


to compute and implement digitally. 


B. WINDOWS USED IN THE SIMULATION 

The windowing schemes that will be implemented in the 
Simulation program DFTWIND.m at Chapter V are described 
below. Figure 31 is the input time domain weighting 
sequence for N = 50 for each of the windows in the 
following. [Ref. 9:p566 to 574] 

1. Rectangular Window 

The rectangular window is one where the inputs are 

uniformly weighted. This is the same as an unweighted 
window as developed in Chapter II and simulated in Chapter 


III. The input weighting sequence is therefore 


w(n) = rect,(n) =n te ees ( 229) 
=0, otherwise. 
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The discrete frequency domain representation of this time 


domain weighting sequence is 


7 a Sb gs 
aaa N sin(2) | a 
N 


where D,(k) 1s known as the Dirichlet kernel. 
2. Bartlett (or Triangular) Window 
This 1S a weighting scheme in which the inputs are 
memoemced by a triangular function. This allows for a 
smoother transition from zero to one for the input weights 


and reduces the sidelobes at the expense of broadening the 


mainlobe. The input weighting sequence is therefore 
win) = (1-218L) rect y(n) . (34) 


The discrete frequency domain representation of this time 


domain weighting sequence is 
_ K))2 
F(k) = [Dn (>) ] . >) 
2 


3. Hanning (Von Hann, or Raised Cosine) Window 
This is a weighting scheme in which the inputs are 
weighted by a raised cosine function such that the magnitude 
of the input at either end of the window is zero. This 
allows for a smoother transition from zero to one than for 


the Bartlett window so the sidelobes are reduced even more 
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However, this is at the expense of broadening the mainlobe 


again. The input weighting sequence is 


PS! 
N 





wih) = racos 


pole 


) Preece, (am (369 


The discrete frequency domain representation of this time 


domain weighting sequence is 
F(k) = Dy(k) += [DeCk-1) DEG (un 


4. Hamming Window 
This is a weighting scheme in which the inputs are 

weighted such that a zero will occur at the peak of the 
first frequency domain sidelobe. This splits the first 
Sidelobe into two smaller sidelobes and is a very good 
method of reducing the appearance of a target report in 
adjacent sidelobes. This gain is paid for with an 
increased processing cost as N needs to be large for the 
best effects. Also, the mainlobe is broadened even further. 
However, the greatly reduced sidelobes may well be worth the 


cost. The input weighting sequence is 


Pel 





w(n) = [0.5435+0.4565 cos{ )] rect =Gae (38) 


The discrete frequency domain representation of this time 


domain weighting sequence is 


F(k) = Dy(k)+0.41996 [D.(k- 1) 290 eom ee (39) 


50 


5. Blackman Window 

The Blackman window is an extension of the Hamming 
menadow in that it manipulates more terms of the Dirichlet 
kernel for greater sidelobe reduction. [Ref. 9:p573] Thus 
the sidelobes are reduced even more than under the Hamming 
window and the mainlobe broadened even more. This gain is 
paid for with an increased processing cost over the Hamming 
as more terms of the Dirichlet kernel must be used. Also, N 
needs to be large for the best effects. Ibs abichelbie 


weighting sequence is 


Wi np) = Ome +0 .5 eos ( 222) 





4mn 


+0.08 cos ( 
N 





Mi Eaeiele (ia): 


The discrete frequency domain representation of this time 


domain weighting sequence is 


F(k) = Dy(kK)+0.52083 [Dy(k-1) +Dy(k+1) 
+0.83333 [(D,y(k-2) +D,(k+2)]. 


C. COMPARISON OF WINDOWS USED IN THE SIMULATION 

The following table is to be used as a rough guide when 
selecting which window function is likely to best suit the 
task. The processing cost is progressive as the windows are 


listed from least to most effective for reducing sidelobe 
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(e) Blackman window 
Figure 31 : Time domain input weighting of the 


various window schemes that are to be simulated in 
Seeeeer V for N = 50. 
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levels. The increase to N required to compensate for the 
mainlobe broadening can be estimated as the mainlobes of the 
windows are given as being a per cent increase over that of 
a mainlobe for a rectangular window using the same N 
samples. Table 1 shows these results for the window schemes 


considered in this thesice 


TABLE 1: CHARACTERISTICS. OF VARIOUS Wl bews 


ei 
LEVEE (ineds) BROADENING 

SS ee 

| wanning | ew 

eee | ee ere 


Blackman 






1S 
ps 
o\o 







Other criteria not given in Table 1 that affect the 


choice of window scheme used are the rate at which the 
sidelobes fall off and the processing cost of the scheme. 
The rate at which the sidelobes fall off for each of the 
window schemes described above can be seen by using the 
Simulation program at Chapter V. The processing cost can be 
estimated by noting how many iterations of the Dirichlet 
kernel need to be calculated for each scheme. Clearly the 
Hamming and Blackman windows require more processing time. 
It 1s worthwhile noting that the Hamming scheme gives both 


improved sidelobe reduction and a reduced mainlobe 
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broadening over the Hanning filter. This may well make the 
greater processing time required for the Hamming window more 
acceptable. 

The simulation of these window schemes can be found at 
Chapter V, and the design process to achieve a suitable 
window scheme with the necessary sidelobe response and an 


acceptable filter width can be addressed. 


55 


Ve SIMULATION OF A WEIGHTED DFT FILTER BANK 


This chapter describes the methodology of the simulation 
of the response of a weighted DFT filter bank to a complex 
input sinusoid, looks at an overview of the program 
structure, describes how to use the program and shows how to 
interpret the results. 

The weighted DFT filter bank simulation program will 
generates a filter bank showing the magnitude response to a 
complex sinusoidal input for all unambiguous frequencies. 
The following filter bank parameters can be varied in the 
Simulation: 

® N = the number of points required in the DFT 
@® PRF = the pulse repetition frequency of the radar 
® k = the filter designator (if required) 


ef... = the doppler frequency of the target (if required). 


A. PURPOSE OF THE SIMULATION PROGRAM 

The equations implemented by the subprograms are fora 
quadrature sampled complex sinusoidal input. The frequency 
domain equations for the magnitude response of the windows 
to be processed were specified in Chapter IV. The increased 
processing time required by the computer for either the 


Hamming or Blackman windows will be evident. 


SiG 


B. PROGRAM STRUCTURE 


Miewoverall program Structure Can be seen in Figure 32. 


The structure is that of a main program (DFTWIND) which 


calls various subprograms (SINGWIN, WINBANK, ... etc) as 


required. There is an error message (ERRMSG) which detects 


an improper data entry, notifies 








DFTWIND.m 





_ SINGWIN.m WINBANK.m 





- ERRMSG.m 


the operator and restarts 


TGTWIN.m 





EXIT 


Figure 32 : Overall program structure for a weighted DFT 


paebeer bank simulation 


the process. It is possible to exit the simulation from 


either the main or subprograms. 

Since it allows maximum operator 
robust format that is data entry 
advantage of this format is that 


to change one of the subprograms 


=e 


This structure was used 
CONGO lnanGeIMpUt.) 9G 1S a 
error tolerant. The other 
should an operator desire 


or add another subprogram 


the system can accommodate this. Any errors made in this 
process will affect only the subprogram in question. The 
rest of the simulation will remain operational while the new 


addition is being debugged. 


C. OPERATING INSTRUCTIONS AND GUIDE 
1. Initial Setup 
To use this simulation the following MATLAB.m files 
must be copied into the user sub directory. Copies of these 


programs are available from Dr GS Gill at USNPGS. 


® DFTWIND.m is the simulation control program. 


@® SINGWIN.m is the simulation subprogram responsible for 
displaying the magnitude response of a single filter 
across all unambiguous frequencies in the filter bank 
for the various windows. 


@® WINBANK.m is the simulation subprogram responsible for 
displaying the full magnitude response for all filters 
in the filter bank for all unambiguous frequencies in 
the filter bank for the various windows. 


® TGTWIN.m is the simulation subprogram responsible for 
calculating and displaying the filter response for a 
target at a given doppler frequency. It can accept 
ambiguous target frequencies and the filter response 
will be for the alias frequency that lies in the filter 
bank. The filter that this alias frequency lies in can 
result in a positive doppler being identified as 
negative or vice versa. 


® ERRMSG.m informs the operator when an entry error has 


been made and restarts DFTBANK after clearing the system 
of all variables. 
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This simulation will work on a DOS 286 machine, but 
a 386 or higher is preferable for increased speed. To start 
the simulation type DFTWIND [ENTER] at the MATLAB prompt. 

2. Printing Graphical Outputs 

The assumption has been made that the PRINT SCREEN 
key will be used to print out plots. To use this key the 
DOS graphics command must have been issued at the DOS prompt 
before any printing is attempted. C:\>> graphics [ENTER] is 
the necessary command. 

Should it be preferable to create meta files for 
each plot the operator has two options. He can either 


® Modify all plot options in the subprograms by adding a 
meta command and file name after the plot command, or 


® Press the Ctrl-C keys to exit the simulation once the 
plot is on the screen. This will return him to the 
MATLAB prompt and he can issue the meta command from 
there. To restart the simulation type DFTWIND [ENTER]. 
3. Program Options 
After the DFTWIND command has been issued to start 
the simulation then the first screen seen by the operator is 
at Figure 33. The second screen the operator will see 1s at 
Figure 34. At this point the operator must choose which 
type of system response that he wants to manipulate. 
Once this has been done the main program transfers control 


to the subprogram selected in Figure 35. The second screen 


is at Figure 36 and the third screen at Figure 37. All of 
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This program will allow you to create weighted 
Doppler Filter Banks for complex sinusoidal inputs. The 
window functions you can use to weight the input data 
are; 


He RECTANGULAR WINDOW 

Pe BARTLETT (or Triangular) WINDOW 

am HANNING (Von Hann or raised cosine) WINDOW 
4 


: HAMMING WINDOW 
ae BLACKMAN WINDOW 
The options are to see the full response of a single 
filter over the full spectrum of non aliasing 
frequencies, to see the full response of a DFT bank of 
filters or to see the magnitude response of each filter 
in a DFT bank to a target occurming at a particule 


frequency. 
TOUCH ANY KEY TO@¢cONT Mise 





Figure 33: First Screen of DFTWIND main control preqiaw 


must now choose which type of response you wish 
to see. 
To see a Single filter response, ENDER 
To see a weighted DFT filter bank, ENTER 2 
To see target values for each filter, ENTER 3 


TO EXIT THIS PROGRAM, ENTER 0Q 
ENTER YOUR CHOICE NOW 
CHOTCE 





Figure 34: Second screen of DFTWIND main control program 
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IF YOU WANT A RECTANGULAR WINDOW, CHOOSE 1 
IF YOU WANT A BARTLETT WINDOW, CHOOSE 2 
IF YOU WANT A HANNING WINDOW, CHOOSE 3 
IF YOU WANT A BLACKMAN WINDOW, CHOOSE 4 
IF YOU WANT A HAMMING WINDOW, Cages k 5 


Then you must decide if you want the response of the 


window chosen to be plotted in comparison with a 
rectangular window or by itself. 
The WINDOW Function you have chosen is # = 


WINDOW = 3 


Compare with a Rectangular Window? Yes = 1, No = 2, 
| compare # = 2 


Figure 35: First screen of a subprogram 





YOU MUST NOW ENTER SOME OF THE NECESSARY PARAMETERS. | 
These parameters are the number ‘N’ of DFT points. 
This is also the number of filters there will be in the 
s filter bank. 
| Next you must choose the PRF you wish to consider. 

Then you must choose the ‘’k’ offset for the DFT to 
determine which filter in the DFT filter bank you will 
select. The choice of ‘k’ is an integer value usually 
meem 0 CO N-1, but you are not limited to this selection 

Finally you must choose whether you wish the display 
of the magnitude response to the DFT to be in normalized 


@enor ratio form. 


BOUCH ANY KEY TO CONTINUE 





Figure 36: Second Screen of a subprogram 


Number of DFT points you wish is, N 
The PRF in Hz you wish is, PRF = 


The Filter K # you wish to use is, 
For a dB plot choose 1, for a ratio plot choose 2, 





Figure 37: Third Screen of a subprogram input request 


Saal 


the screens for the subprograms are very Similar, the only 
difference being whether" "k" or “EU Gre” pequser = cess 

First the choice of which window scheme to examine 
must be made. Then the decision whether or not to plot the 
selected window in comparison with a rectangular window or 
by itself must be made. The plots in the simulation are 
done in different colours and line styles so they are easily 
distinguished on the screen. Green is for the rectangular 
window and red for the window being studied. The colours 
will not show up on the printout of the plot so it can be 
difficult to differentiate between line styles. The chou. 
of a dB or ratio plot on the third screen of a subprogram 
determines whether the magnitude scale on the y axis will be 
plotted for ratio values or in a decibel format. When 
attempting to determine the sidelobe level response the dB 
format is preferable. For a more general overview the ratio 
format is better. 

Once all of the graphical outputs and results have 
been shown, each subprogram will terminate with the screen 
at Figure 38 where the operator can choose to do the same 
subprogram for the same type of input to the DFT (the 
SINUSOID variable remains constant), return to the main 
program DFTBANK to try a different simulation or exit the 


Simba ome 
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Should any data entry error have occurred that the 
inherent MATLAB error routines did not recognize, than the 
Simulation will detect it and the error message at Figure 39 
will be displayed. All variables will be cleared and the 


main program DFTBANK will be restarted. 


— — ——————e EE ———— 


meeevUol NOW CHOOSE TO EITHER BO THIS CHOICE AGAIN, 
Par TO THE MAIN MENU OR EXIT THE SIMULATION 







DO THE SINGLE LOBE AGAIN, CHOOSE. 
TO RETURN TO THE MAIN MENU, GHOOSE 2 
mo EXTT THE SIMULATION, CHOOSE 0 


Figure 38: Final Screen of a subprogram 


YOUR CHOICE DOES NOT CORRESPOND TO THOSE AVAILABLE. 
THE PROGRAM WILL BE RESTARTED FOR YOU. 
REMEMBER YOU MUST CHOOSE OQ TO EXIT. 


PRESS ENTER TO CONTINUE... 





Figure 39: Error message subprogram screen 


D. INTERPRETATION OF RESULTS 

The output of the simulation is mostly graphical in 
nature Since graphical displays are readily assimilated and 
provide a more intuitive understanding of the process. 
Results from each of the subprograms will be presented. 
These plots will be analyzed to demonstrate some of the 
Capabilities of the simulation program and to provide worked 


examples as a reference guide. In each case the required 
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data inputs will be provided and the selection process is 
described. 
1. Single Filter Response 
The SINGWIN subprogram displays the magnitude 
response of a single DFT filter weighted according to choice 
across all unambiguous frequencies in the filter bank. 
These unambiguous frequencies are for |f |s PRF/2. While it 
is the frequency range of the mainlobe magnitude response of 
the filter that defines the frequency coverage of the 
filter, however there can be a response to any frequency, 
even from the sidelobes of the filter. If the response for 
a frequency were from the sidelobes of the filter and it 
exceeded the detection threshold set for that filter, a 
false target report for that filter would be given. 
The operator can vary the PRF and N in the 

Simulation to achieve a filter response with an acceptable 

@e filter width, 

@® frequency resolution, and 

® doppler frequency coverage. 
The data at Table 2 is good for a quick comparison of the 
relative windows and can help in determining the most likely 
windowing scheme to satisfy a radar design requirement. 
Table 2 shows these results for the windowing schemes that 


were addressed in this ™tnesic. 


64 


MAINLOBE 
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Blackman | 89 





In the following example the inputs to the 


Peemulation were 


® CHOICE = 1 for a single filter response 

@® WINDOW = 2 for a Bartlett window 

® Compare = 2 so it will not be compared toa 
rectangular window 

@®N = 8 

® PRF = 2000 Hz 


@®k = 0 for the F(k) magnitude response 
Smaeeratio = 2 for a ratio plot 
The plot resulting from these simulation choices is 
at Figure 27. Had neo = 1 been selected than the dB 
ratio magnitude response at Figure 28 would have been 
plotted. 
2. Weighted DFT Filter Response 
The WINBANK subprogram displays the mainlobe 
magnitude response of all filters in the filter bank for all 


unambiguous frequencies. The choices are the same as for 
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SINGLOBE, with the exception of "k" since the mainlobe 
magnitude responses for k=0 to N-1, k ¢€ I are all plotted 
automatically. 

In the following example the inputs to the 


Simulation were 


I 


® CHOICE = 2 for a full filter bank response 


| 


@® WINDOW = 4 for a Hamming window 


® Compare = 1 so it will be plotted versus a 
rectangular window 


@® PRF = 2000 Hz 
® dB ratio, =) itor -aede- peer 
The plot resulting from these simulation choices is 
at Figure 41. From it the filter width at crossover of the 
mainlobes, straddle loss and other performance measurements 
in comparison to an equivalent rectangular window DFT can be 
noted. This demonstrates where the detection threshold can 
be set for the weighted window in comparison to the 
unweighted rectangular window if it is desired to 
® avoid declaring targets in more than one mainlobe of a 
filter, then the threshold must be set higher than the 
mainlobe crossover point, or 
@® avoid declaring targets from the sidelobe of a filter, 
then the threshold must be set higher than the highest 
sidelobe. 


This gives the operator a good feel for the change in 


dynamic range for the different windows over which he can 
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Figure 40 : Magnitude response for the F(0) filter of a 
Bartlett weighted DFT 
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Figure 41 : Magnitude response for the full filter bank using 
a Hamming weighted DFT in comparison to an unweighted 
Beoeramgular DFT 
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detect targets that vary in signal strength and not affect 
his false alarm rate. 

The graph is difficult to analyze when it is %pe@mee 
in black and white, especially if the additional information 
from the rectangular window for reference purposes is on it 
as well. The graph is much more intelligible when viewed on 
screen as the different colours make interpretation easier. 

3. Filter Bank Response to a Specific Target Frequency 

The TGTWIN subprogram is the same as the WINBANK 
Subprogram, except that it will require an additional input 
for the doppler frequency the target or clutter is to appear 
at and that all plots are automatically of a dB nature. If 
the doppler frequency given does not satisfy the Nyquist 
criteria then the alias frequency (f,) will be processed by 
the filter bank. The filter that the alias frequency lies 
in has no direct relation as to whether the original doppler 
waS a positive or negative doppler frequency. Thus, this 
can result in a positive doppler being identified as 
negative or vice versa. 

In the following example the inputs to the 
Simulation were 


@® CHOICE 


3 for a target repore, 


@® WINDOW 


4 for a Hamming window, 


® Compare = 1 so it will be plotted versus a 
rectangular window, 


® N = 16, 
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@® PRF = 2000 Hz, 
Seaemeatio = 1 for a dB plot, and 
See, - 220 Hz. 


The plot resulting from these simulation choices is at 


Figure 30. The sidelobe improvement of the Hamming window 


TARGET RETURN FOR f-—tgt = 356 We with co PREF @w» 2000 Hz and N = 16 





Normalized Magnitude 





Filter Bonk # 


Figure 42 : Target report for a Hamming window (0) in 
comparison to that from a Rectangular window (x). The 
Hamming window shows the steep sidelobe drop off that 


makes it useful 


over that of an unweighted DFT is evident, though it is 
worth noting that the sidelobes adjacent to the mainlobe do 


not see the same improvement that the other sidelobes do. 
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fe PULSE COMPRESSION TECHNIQUES 


Pulse compression (PC) is a technique in which a long 
pulse that has been modulated in either frequency, amplitude 
or phase is transmitted. The received signal is then 
processed ina matched filter to compress it to a narrow 
pulse. Doing so has the following advantages and 
disadvantages [Ref. 6:p420 to 434], [Ref. 4:p123 to 156] and 
FREE olin. 


® A long pulse can be used to increase the average 
transmitted power. 


® The short pulse in use after compression is desired for 
range resolution, accurate range measurements and 
discrimination of multiple targets. 


@® The radar is less vulnerable to ECM that cannot generate 
the same coherent waveforms as the radar. 


® The target signal experiences a processing gain while 
the noise in the system does not. 


@® Some pulse compression waveforms are sensitive to 
doppler effects and this may degrade the overall 
performance of pulse compression. 

@® The minimum range at which the target can be detected is 
increased for a radar with pulse compression. 

A. RADAR RANGE AND RESOLUTION 
Two of the key characteristics of a radar system are a 
long detection range and a narrow range resolution (AR). 


The requirement for a long detection range is self evident. 


The reason to have a narrow range resolution is to resolve 


1G 


two targets that are close together. Both of these 


characteristics are linked by a dependence on the pulse 


length of the transmitted waveform. 


1. Range > 
The range dependence on the pulse width (T) can be 
seen from the radar range equation as follows: 
ae DeeeG*A*o | wy 
(re) ei Bia 
nor 
Pave = a = (43) 

where 

® R = the radar detection range, 

@® Prayg = the average power of the transmitted pulse, 

® G = the gain of the monostatic radar’s antennae, 

@® \ = the wavelength of the transmitted waveform, 

@® o = the Radar Cross Section (RCS) of the target, 

® F, = noise figure, 

® k = Boltzman’s constant, 

® T = the temperature the radar is operating at, 

@® B = the filter bandwidth, 

e L. = the system losses in the radar, 

® Pohax = the peak power the radar is operating at, 

® T = the pulse width of the transmitted waveform, and 


a 


@® PRI = the pulse repetition inven 


Thus if all other factors in Equation 42 arewiiess 
constant except for P,y,, the only way to increase the 
detection range is by increasing either P,.,, or T. There 
are practical limits on the peak power a radar transmitter 
can generate. With all other factors in Equation 35 
remaining constant, the only way to effect a further 
increase in the range of the radar is by increasing T. 

2. Range Resolution 
A good range resolution (AR) is desirable for the 


following reasons: [Rete 6: p42)) 


@® Better range resolution leads to better range accuracy. 


® When attempting to separate multiple targets it is 
necessary to separate them in either range, azimuth or 
elevation. Given the 3 dB beamwidth of antennas 
Suitable for search radar applications a good angle 
resolution in either azimuth or elevation is difficult 


to achieve. Thus it is preferable to separate multiple 


targets in range rather than in angle. Therefore the 


smaller the range resolution is, the better the ability 


to separate multiple targets in range will be. 


® The amount of clutter a target has to compete against 
varies directly with AR. Thus the smaller AR is, the 


less clutter there is to raise the CFAR detection level. 


Therefore it will be possible to detect smaller 
magnitude targets. 


The equation for the calculation of a radar's AR is 


as follows: 


ee 


AS ee (44) 


Clearly AR depends only on T, hence a shorter T is 
desirable. Therefore there is a conflict between the 
requirements for a high range resolution (small T) and a 
large detection range (large T). One way to try to meet 


both of these requirements is through pulse compression. 


B. PULSE COMPRESSION 

Pulse compression (PC) is a scheme to resolve the 
conflicting requirements of narrow AR and long detection 
range. In pulse compression a long pulse (T) composed of 
several sub pulses (7T) which have been modulated in phase, 
frequency or amplitude is transmitted. The received signal 
is correlated with the replica of the transmitted pulse. 
This results in the correlation having a high impulse like 
output when there is a return from a target and a low output 
in the absence of a target. The long transmitted pulse, 
Memersting Of all the sub-pulses, increases P,,, to its 
Seeeeeed value. Thus P,,, continues to be dependent on T. 
However, the narrow pulse output of the correlator will 
satisfies the AR requirement and is not dependent on T. 

1. Pulse Compression Ratio 

The relationship amoung T, 7 and the PRI is seen in 

Figure 43. If the transmitted pulse width is T and the 


subpulse width is 7, then, 


Ue 


4 
TL = NG. se 


The pulse compression ratio (PCR) is 


L (46) 


PC = 
_ 


The pulse compression achieved in a real world system is 
usually some value less than or equal to the PCR. Any 


losses due to doppler or other effects are referred to as 


the pulse compression loss factor. 





Magnitude 
-—~- 4 fi —~-—— 
Tt | v 
5 | 
PRI Time 
(T). “wie 


Relationship between the pulse width 


Figure 43 
and the PRI for pulse compression 


subpulse width (rT) 


2. Range Resolution 
The range resolution under pulse compression 1s a 


function of the subpulse length. Therefore, for pulse 


compression coding schemes 
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Now Piyg « T while AR « 7; thus with pulse 
compression, it 1s possible to achieve both a good detection 
range and a good range resolution. 

3. Pulse Compression Processing Gain 

Processing gain is the amount by which the SNR at 
the output of PC improves with reference to the SNR at the 
input. Noise at the input 1s random and has a low 
correlation with the replica of the transmitted pulse. 
However, the target return has a high correlation and 
receives a processing gain equal to the achieved pulse 
compression. Therefore, the SNR is improved by a factor 
equal to the pulse compression. This is one of the 
Significant advantages of using pulse compression 
techniques. 

4. Eclipsing Loss 

When a radar pulse is being transmitted, it 1s not 
possible for the receiver to operate. Any target information 
that appears during this time interval is lost. For a radar 
uSing pulse compression, the length of the transmitted pulse 
1S increased and more eclipsing loss will result. This will 
affect the minimum range (R,,,) at which the radar can detect 


targets. Also, if the radar is designed to detect targets 


TD 


beyond the unambiguous range (R,), the eclipsing loss is 


repeated at each multiple of R,. 


Cz PULSE COMPRESSION PERFORMANCE CRITERION 
The measures of PC performance used in the simulation 
programs at Chapter VII are defined as follows; 
1. Peak to Sidelobe Ratio 
The peak to sidelobe ratio (PSL) is the ratioue@eene 
highest sidelobe to the mainlobe. High sidelobes tend to 
cause false alarms when they cross the detection threshold 
Setting. To reduce the false alarm rate, the detection 
threshold must then be raised, thereby desensitizing the 
radar. It can be seen that a small target could thus be 
missed by the radar. 
2. Integrated Sidelobe Ratio 
The integrated sidelobe ratio (ISL) references the 
effect of the sum of all the sidelobes to the mainlobe ae 
is desirable to have the maximum energy possible in the 
mainlobe, where it will be used for detection purposes since 
any energy in the Beene. is discarded insofar as target 
detection is concerned. 
3. Pulse Compression Ratio 
The pulse compression ratio (PCR) is the ratio of 


the transmitted pulse width to the compressed pulse width. 
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4. Doppler Sensitivity 

When the transmitted radar pulse encounters a moving 
target, a doppler shift occurs in the reflected pulse. A 
doppler shifted waveform from a moving target may not 
compress to its full potential and thus provides only a 
limited processing gain. This effect varies from code to 
code. Codes which can continue to compress effectively in 
spite of doppler shifts are known as doppler tolerant codes 
and thus are preferable. 


Target motion causes a shift in frequency given by; 





f, = | (48) 


The corresponding change in phase is given by: 


Ore 2a ets FOL n- 1 -N, ne ©, (49) 
where 
® V, = the relative doppler velocity between the target 
and the radar, 
® f, = the doppler frequency resulting from V,, 
® \ = the wavelength of the transmitted pulse, 
® ¢, = the doppler phase shift resulting from Vj, 
®n = the subpulse number, and 
® 7 = the subpulse width. 


Since the doppler phase shift is progressive from 
subpulse to subpulse, this may pose a serious problem to 


phase coded pulse compression schemes. It can cause 
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reductions in the processing gain and increases in the 
Sidelobes. If a code is not doppler tolerant then there are 
schemes available to attempt to reduce the influence of the 
doppler effect. One of these schemes is described in the 


Simulation of pulse compression in Chapter VII. 


D. TYPES OF PULSE COMPRESSION 
The three basic types of modulation schemes in use are 
amplitude, frequency and phase modulation. 
1. Amplitude Modulation 
This modulation scheme is almost never used since 
the amplitude modulation of a pulse train is so easily 
distorted by numerous of external factors. [Ref. 6:p420] 
2. Frequency Modulation 
Historically the most common pulse compression 
schemes used have been those of the frequency modulated 
form. In this scheme the frequency modulation of the 
carrier is varied throughout the overall pulse length. This 
variation can be either linear or non-linear and it can be 
implemented in either a continuous or discrete form. 
Therefore the four subsets of frequency modulation are; 
(Ref. 3:p10.4], [Ref. 2:p9 to 25) continuveuSswiiaea 
frequency modulation, discrete linear frequency modulation, 
continuous non linear frequency modulation and discrete non 


linear frequency modulation. 
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a. Continuous Linear Frequency Modulation 

The most common type of this form of modulation 
as linear frequency modulation (LFM). In LFM, there are no 
actual subpulses since the waveform is continuous. The 
frequency of each transmitted pulse of duration T is 
increased or decreased at a constant rate over the pulse 
length as shown in Figure 44. Tne Mat@eed filter for this 
waveform has a "lagtime versus frequency" characteristic 
that is the opposite of the transmitted pulse. Thus when 
the returned signal is passed through the matched filter, 


Bmewcrailing portions of the signal "catch up" to the 


Frequeocy LINEAR FREQUENCY MODULATION 
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Figure 44 : Linear frequency modulation scheme 


leading edge, and a compressed output with a greater 
amplitude and smaller pulse width results. [Ref. 1:p217] 
The amplitude gain over the transmitted signal is the PCR, 


and the compressed pulse width corresponds to Tt=T/PCR. The 


Us: 


design specifications for LFM are as follows; [Ref. 2:p13], 
(Ref. 12] 


@e Af is the minimum resolvable frequency difference, and 
it is defined as being Af=r". 


@® In order to achieve Af the modulation rate (F) Gveremeae 
pulse length must be at a rate no less than F=Af/T. 


@® Clearly for a fixed T, the width of the compressed pulse 
1s determined by the modulation rate. The greater the 
rate of change for f as seen in Figure 45, then the 
smaller +r will be. 

@® The overall pulse is formed by starting with a pulse 
equal to 7. Then the frequency is linearly changed at a 
rate where T=7r+kAf. Here k is the dispersion 
coefficient in units of sec/Hz/sec from the dispersive 


delay line that forms the pulse T. Since in almost all 
cases kAf » 7, the PCR is considered to be 


pcr. = + = 21802 LENE. (50) 
4 & a % 


Target doppler causes an error in the range 
measurement for the LFM waveform. However, there is 
no other significant adverse effect due to the doppler 
shift. Range errors can be corrected if the target doppler 
shift is known. 

The other problem of LFM is that the time 
Sidelobes of the compressed waveform are large enough to 
affect radar sensitivity and the constant false alarm rate 
(CFAR). Large targets have significantly large sidelobes 
which interfere with the detection of small targets. If the 
detection threshold is set to detect small targets, then the 


Sidelobes of large targets are reported as targets as well. 
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ftmchne thresholdwais raisedwhagh enough to avoid this 
problem, than small targets go undetected. 


The chief advantages of the LFM waveform are its 


doppler tolerance and the ability to achieve very large 
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Figure 45 : Linear frequency modulation scheme undergoing a 

doppler frequency shift 

PCRs. The time sidelobe problem can be reduced through the 


use of weighting schemes, such as those introduced to reduce 
the frequency sidelobes of the DFT filter bank in Chapter 
Ty. 
b. Discrete Linear Frequency Modulation 
In discrete linear frequency modulation the 
frequency is constant over a given subpulse interval 7 and 
then instantaneously changed by an amount Af for the next 


subpulse, unlike the continuous linear change for LFM as 


ga 


seen in Figure 46. This is done for a total of N subpulses. 
The best frequency step interval found has been to have 
Af=7' which avoids undesirable pulse compression waveform 


sidelobes. Digital LFM has the same characteristics as LFM. 


| Rein 2:pisi 


Frequency 


Time 
Pulse Width 
Figure 46 : Digital linear frequency modulation scheme 
c. Continuous Nonlinear Frequency Modulation 

Nonlinear FM waveforms can be used in modern 
radar systems. They have some advantages over linear FM 
waveforms, the most important being their greatly reduced 
time sidelobes after compression. However, their chief 


disadvantage of expensive hardware implementation outweighs 


this gain, except for very large frequency invariant 
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systems. This waveform is used in FAA radars for enroute 
commercial aircraft surveillance. 
d. Discrete Nonlinear Frequency Modulation 
Digital waveforms are easier to implement and 
process in modern computers. However, it is doppler 
intolerant and just as expensive to implement as the 
continuous nonlinear FM. 
3. Phase Modulation 
This 1s another form of pulse compression modulation 
in practical use. Its advantages include better time 
sidelobes than for FM, the ease of generating the waveform 
and use of the same hardware to generate different length 
codes. There are two basic types of phase modulation in use 
today, and they are binary phase or polyphase modulation. 
Bees 10.15 to 10.27), [Ref. 2:p17], (Ref. 12) 
a. Binary Phase Modulation 
In binary phase modulation, the phase changes to 
the waveform are in either 0 or mw increments from subpulse 
to subpulse for a total of N subpulses. Since T=N7, then by 


definition the PCR is 
PCR = = =ONE (51) 


The time sidelobes from the correlator for binary codes need 
to be minimized as much as possible, just as in the FM case. 


Thus, based on this criterion, the binary phase modulation 
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schemes are broken into two subclasses, namely Barker and 
pseudo random codes. 

(1) Barker Codes. The Barker codes are known as 
perfect codes because their time sidelobes are all equal and 
of unit magnitude. Thus the waveform has the minimum 
possible energy contained in the ISL, and all sidelobes are 
at the same PSL. [Ref. 2:p17] However, they are only known 
to exist up to a maximum length of 13. After this length, 
it 1s possible to use combined Barker codes for a greater 
PCR in which the Barker subpulses are themselves coded 
within another Barker code. While this extends the PCR, the 
PSL remains the same as that of the lowest order Barker 
sequence used in the combined coding, and the ISL worsens. 
The main problem with Barker and combined Barker codes is 
their low tolerance to doppler shifts as the performance of 
the received waveform in the matched filter 1s seriously 
affected; PCR, PSL and ISL suffer. If the transmitted pulse 
length is short, the adverse doppler effects are kept within 
tolerable limits for the Barker codes. However, when the 
combined Barker is used to increase the PCR, its greater 
length leaves it vulnerable to the adverse doppler shift 
effects. The Barker codes are not particularly immune to 
the effects of ECM. Once the waveform is detected and 
identified it will be susceptible to ECM that has the 


capability to generate coherent phase coded waveforms. 
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(2) Pseudo Random Codes. If the binary sequence 
transmitted is a purely random sequence of 180° phase shifts 
for N subpulses, then the PCR is still be equal to N. 
However, the PSL and ISL characteristics of the waveform are 
most likely poor due to high time sidelobes resulting from 
the correlation process. For this reason, truly random 
codes are seldom used. However, an attempt to approximate 
random phase codes whilst achieving the best ISL and PSL 
possible for that length of code is made because of the high 
ESM tolerance of random codes. The most common form of the 
pseudo random codes is the one known as maximal length 
sequence, where an "n" stage shift register is used with 
Feedback to generate the code. The code will be of length 
N = 2°-1, and the code created by this process depends on 
the "n" initial conditions set in the shift registers and 
the form of the feedback loop. The feedback loop can be a 
modulo adder, XOR gate or any other logical device. Figure 
47 illustrates this process for a modulo 2 adder with inputs 
from registers three and four for a four stage register. 
Bees: p10.18 to 10.20], (Ref. 6:p428 to 431] 

Clearly a great number of codes can be 
generated for any given N. Many of these codes have poor 
ISL and PSL performance on compression and are discarded. 
The code or codes which for any given N show the best ISL 


and PSL performance are known as the optimal codes. A great 
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many optimal codes are known to exist and are in use today. 
For example, for N=40, there are 114 known optimal codes, 
and thus the choices for ESM suppression are widespread. 
[Ret 3olleuais | 

The same problems that occur with binary 
codes in respect of doppler phase tolerance also occur with 
pseudo random codes. In addition, since the time sidelobes 
of pseudo random codes are not of unitary height the doppler 
performance of pseudo random codes is usually worse than for 


the Barker codes. 


Modulo 2 —————- 
Adder |__ 


Pigure 47 : A pseudo random code generator using an adder with 
four initial conditions to be set. A maximal length 15 bit 
pseudo random code can be generated. 


b. Polyphase Modulation 
Rather than having the subpulses use phases of 0 


or 180°, the various polyphase schemes use multiple values 
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between 0 and 360° for the subpulse phase modulation. The 
schemes for the creation of the various polyphase schemes 
vary Significantly in scope, characteristics and results. 
Those polyphase codes based on nonlinear waveform generation 
have the same problems as do the nonlinear FM codings 
discussed earlier. For this thesis, only some of the 
polyphase schemes based on the LFM modulation scheme as 
adapted for polyphase implementation are examined. These 
codes are the Frank polyphase code and the P3 polyphase 
code. Since both of these codes are derived from the LFM 
waveform, they possess all of the properties, limitations 
and advantages of LFM modulation. Since they can be 
implemented digitally, they are easy to use in modern 


meemeae (Ref. 3:p10.25], (Ref. 2:p15] 
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i gn SIMULATION OF PULSE COMPRESSION 

This chapter describes the methodology of the simulation 
of pulse compression for a variety of modulation schemes and 
examines one of the methods of correcting for doppler shifts 
in pulse compression. It also looks at an overview of the 
program structure, describes how to use the program and 
shows how to interpret the results. 

The pulse compression simulation program will generate 
the output of the correlator for the chosen modulation 
scheme taking into account any doppler effects present. The 
shape of the correlation is displayed graphically, along 
with the PCR, PSL, ISL and the pulse compression achieved by 
the correlator. The following parameters can be varied in 
the simulation: 

® The type of modulation scheme to be used, 
@® N = the number of subpulses in the modulation scheme, 
® PRF = the pulse repetition frequency (if required), 


® RF = the frequency the radar is transmitting at (if 
required), 


® 7 = tau, the subpulse width, 


® V,, = the actual speed of the target (1f required) iam 
® Vi. = the assumed maximum speed of the target (if 
required). 


The modulation schemes that are simulated are: 


® Barker codes from binary phase modulation, 
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the 


Combined Barker codes from binary phase modulation, 


Pseudo Random codes with a modulo 2 XOR gate from binary 
phase modulation, 


P3 codes from polyphase modulation (which has LFM 
characteristics), and 


Frank codes from polyphase modulation (which has LFM 
characteristics). 


By controlling these parameters the student can address 


realities of some of the pulse compression schemes in 


current use. Some of the situations that can be addressed 


ale: 


ake 


The various performance criteria of a given modulation 
scheme can be assessed. 


The performance criteria of a given modulation scheme 
with a target doppler shift can be assessed. The 
doppler shift can either be for a given target velocity 
Or an assumed doppler shift across the specified pulse 
woath . 


The effects that the transmitted radar waveform pulse 
width will have on the eclipsing loss for a specific 
radar can be assessed. This is presented graphically by 
displaying both T and the PRI over three pulses. The 
minimum range (R,,,) that the radar can detect at and its 
unambiguous range (R,) will be given. 

There are numerous other possibilities left to the 
user's initiative. 

PURPOSE OF THE SIMULATION PROGRAM 

The modulation schemes implemented by the sub programs 
detailed below. 

1. Barker Code 


The Barker codes are the only codes which have 


sidelobes of unit magnitude and thus yield the best results 
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for their lengths. They are formed by changing the phagevon 
the transmitted frequency by a phase shift of 0 or mw for 
each 7 interval. These 0 and mw phase shifts are normally 
expressed by taking their cosines and thus representing them 


as binary codes of +1. There are N subpulses and the PCR is 


PCR = = N. (52) 


Although this code is doppler intolerant its short 
length may help shield it from the worst effects. The known 
Barker codes are listed in Table 3. An example of the 
output of the correlator for a Barker 13 code witheus 
doppler effects is at Figure 48. 

2. Combined Barker Code 

This 1s where a Barker code is coded within another 
Barker code to increase the length of the code, and thus the 
PCR. PSL remains the same as that of the lowest Barker code 
used, and the ISL usually worsens. This code is not doppler 
tolerant, and since it is longer than the Barker codes it is 
more susceptible to doppler. It is usually referred to asa 
‘m’ in ‘n’ combined Barker code, where m is the inner code 
and n is the outer code. Thus there are a total of N=m*n 


subpulses, and the PCR is 
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TABLE 32: ~ THE. KNOWN BARKER CODES 
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Beomre 48 : Output of the correlator for a Barker 13 code 
without doppler effects 
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PCR = m+n ="0e (53) 


ae example of a 13 in 13 combined Barker code is at Figure 
3. Pseudo Random Codes 

The pseudo random codes can be generated fora 
length n shift register using a modulo 2 XOR gate. Any two 
of the shift register cells can be selected to be fed to the 
XOR gate. The shift registers to be used are referred to as 
PN(3,4) when the 3" and 4" shift register cells are to be 
used as feedback to the XOR gate. The number of "n" initial 
conditions must also be specified. The maximum length for a 
unique code is N=2°-1. This code is useful for ECM 
purposes, but has a poor response to uncorrected doppler 
conditions. Also, it has large sidelobes when compared to 
the Barker codes, which is a major disadvantage. The 
correlation of a n=4, PN(3,4) code with all au1e7e 
conditions set to 1 is at Pugure Soe 

4. P3 Polyphase Code 

There are many variations of this type of code, but 
the one considered here is the P3 code derived from LFM 
(which makes it doppler tolerant). The mathematical form 
for the generation of the code phases for the P3 polyphase 


code is 
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Figure 49 : Output of the correlator for a combined Barker 
ian 13 code without doppler effects 


Pseudo Roaondam Code for n = 4, shifter at p = 3B andaq-= 4 
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Figure 50: Output of the correlator for a pseudo random code 
fMeeweut GOppler effects for 4 initial conditions (all equal 
to 1)and an XOR gate feedback from bins 3,4 
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(7-1)* | 2721S Naeem (54) 


where N is the number of subpulses desired (also the PCR). 
An example of the output of the correlator for the P3 code 
1S At Figunmers i 
5. Frank Code 
The mathematical form for the generation of Frank 


code phases is 
(as) (55) 
N 


where i=j=n, the Frank code number, ne I, and N = n° is the 
number of subpulses present. 

The code is formed in much the same fashion as the 
combined Barker where 1 corresponds to the inner code and j 
to the outer code. The resulting ¢ matrix is reshaped into 
a vector subpulse train. An example of the output of the 
correlator for a Frank code with n=4 and N=16 is at Figure 


52% 


B. PROGRAM STRUCTURE 

The overall program structure can be seen in Figureweas 
The main program (PCRSIM) which calls one of three 
Subprograms NODOPP, DOPP or XDOPP depending on whether 


processing is done without the effects of doppler, with 
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mra@ure S51: Output of the correlator for a P3 polyphase code 
Meee = 169 without doppler effects 
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meee 52: Output of the correlator for a Frank code of n=4, 
Mees N=16) without doppler effects 
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uncorrected doppler or with corrected doppler. There islam 
error message (ERR2MSG) which detects improper data entries, 
notifies the operator and restarts the process. Then the 
Choice of modulation scheme is made, and this can be seen in 
Figure 54, which is an expanded portion of Figure 53. Due 
to the peculiarities of MATLAB, it is necessary to have 
three separate sub programs for each modulation scheme to 
account for the presence of no doppler, uncorrected d@maeees 
or corrected doppler. It is possible to exit the simula 
from either the main or sub programs. This structure was 
used since it allowed the best user interface and is a 
robust format, that 1s data entry error toleramite The other 
advantage of this format is that should an operator desire 
to change one of the sub programs or add another sub program 
the system can accommodate this. Any errors made in this 
process will affect only the subprogram in question. The 
rest of the simulation will remain operational while the new 


addition is being debugged. 


C. DOPPLER CORRECTION SCHEME FOR PULSE COMPRESSION 

When the transmitted radar pulse encounters a moving 
target, a doppler shift occurs in the reflected pulse eae 
target's relative velocity to the radar causes a shift in 


frequency given by: 
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Figure 53 : Overall program structure for PCR simulation 
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Figure 54 : Sub program details of PCR simulation 
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The corresponding change in phase is given by: 


$6. = 20£, epee =" eecans (S79) 


where 


@® V, = the relative doppler velocity between the target 
and the radar, 


@ 

rh 
ro 

it 


the doppler frequency resulting weve 
@e \ = the wavelength of the transmitted pulse, 
@® ¢, = the doppler phase shift resulting from V,, 


@®n the subpulse number, and 


e@ 
~ 
lt 


the subpulse width. 


If the radar platform is itself in motion this cause an 
additional doppler shift as well. However, since the speed 
of the radar platform is known, it can be readily countered 
in the RSP calculations and is therefore not included in the 
following equations. One of the simplest schemes in common 
use 1s to attempt to counter the effect of the doppler shift 
by introducing another shift of equal Magnitude seme 
opposite in sign. Thus if a doppler shift 1s negative, a 
positive shift of the same magnitude is used to cancel it 


Oui: 
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This scheme is demonstrated in the XDOPP.m subprogram 
option of the PCRSIM simulation program. The concept behind 
this correction scheme is as follows; 


@® If the radar itself is in motion this speed will be 
known, and it can be accounted for in the assumed 
Beuneer doppler equation. 


® The highest possible speed (Vy,,)of the target to be 
detected 1s Known or can be assumed. One half of this 
speed is then used in the counter doppler equation at 
Equation 74 for V,. The assumption of one half maximum 
is made since in most situations this will approximate 
the likely speeds encountered. 


@® A decision as to whether the radar will be encountering 
@eening or closing targets 1S critical. If the wrong 
decision is made then the overall doppler shift 
encountered will be worse than if no correction were 
made. If the target’s doppler shift is unambiguous to 
the radar’s filter bank the correct decision will be 
made. However, if the dop pler frequency is ambiguous 
than it is quite possible that the wrong decision will 
be made. 


@® If the correct decision as to an opening or closing 
target has been made this simple technique is 


surprisingly effective in removing the worst effects of 
doppler. 


The equations for this correction scheme are as follows: 


SESH, 


Sy MAK ! 
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which results in the overall doppler shift the waveform sees 
being (if the correct decision as to an opening or closing 


target is made) 


oo 


, 2 2 Wave) oa 


or (1£ the wrong decisionm@is made, 


2 {| Vee 
fy = 2a (60 


The corresponding change in phase uses the new f, from 


Equation 59 or 60 and is given by: 
6, = 2nf, nt fOr m=1-\eaeoee (61) 


where V, 1s the counter velocity and Vy,, is the maximum 


velocity of the target. 


D. OPERATING INSTRUCTIONS AND GUIDE 
1. Initial Setup 
To use this simulation, the following MATLAB.m files 
must be copied into the user sub directory. Copies of 
these programs are available from Dr GS Gill at USNPGS. 
® PCRSIM.m is the simulation control program. 
® NODOPP.m is the simulation subprogram responsible for 
processing a modulation scheme without the presence of 


doppler. 


@® DOPP.m is the simulation subprogram responsible for 
processing a modulation scheme with doppler. 


@® XDOPP.m is the simulation subprogram responsible for 
processing a modulation scheme with corrected doppler. 


@® BARKER1.m is the simulation subprogram responsible for 


processing the Barker modulation schemes without the 
presence of doppler. 
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BARKER2.m 1s the simulation subprogram responsible for 
processing the Barker modulation schemes with doppler. 


BARKER3.m is the simulation subprogram responsible for 
processing the Barker modulation schemes with corrected 
doppler. 


COMBO1.m is the simulation subprogram responsible for 
processing the combined Barker modulation schemes 
without the presence of doppler. 


COMBO2.m is the simulation subprogram responsible for 
processing the combined Barker modulation schemes with 
doppler. 


COMBO3.m is the simulation subprogram responsible for 
processing the combined Barker modulation schemes with 
corrected doppler. 


RANDOM1.m is the simulation subprogram responsible for 
processing the pseudo random modulation schemes without 
the presence of doppler. 


RANDOM2.m is the simulation subprogram responsible for 
processing the pseudo random modulation schemes with 
doppler. 


RANDOM3.m is the simulation subprogram responsible for 
processing the pseudo random modulation schemes with 
corrected doppler. 


POLYPH1.m is the simulation subprogram responsible for 
processing the P3 modulation schemes without the 
presence of doppler. 


POLYPH2.m is the simulation subprogram responsible for 
processing the P3 modulation schemes with doppler. 


POLYPH3.m is the simulation subprogram responsible for 
processing the P3 modulation schemes with corrected 
doppler. 


FRANK1.m is the simulation subprogram responsible for 
processing the Barker modulation schemes without the 
presence of doppler. 


FRANK2.m is the simulation subprogram responsible for 
processing the Barker modulation schemes with doppler. 
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@® FRANK3.m is the simulation subprogram responsible for 
processing the Barker modulation schemes with corrected 
doppler. 

@® ECLIPS.m is the simulation subprogram responsible for 
displaying the eclipsing loss parameters after one the 
modulation subprograms has been calculated. 

® ERR2MSG.m informs the operator when an entry error has 
been made and restarts DFTBANK after clearing the system 
of all variables. 

This Simulation will work on a DOS 286 machine, but 
a 386 or higher is preferable for increased speed. To start 
the simulation type PCRSIM [ENTER] at the MATLAB prompt. 

2. Printing Graphical Outputs 

The assumption has been made that the PRINT SCREEN 
key will be used to print out plots. To use this key the 
DOS graphics command must have been issued at the DOS prompt 
before any printing is attempted. C:\>> graphics [ENTER] is 
the necessary command. 

Should it be preferable to create meta files for 


each plot the operator has two options. He can either 


® Modify all plot options in the subprograms by adding a 
meta command and file name after the plot command, or 


® Press the Ctrl-C keys to exit the simulation once the 
plot is on the screen. This will return him to the 
MATLAB prompt and he can issue the meta command from 
there. To restart the simulation type PCRSIM [ENTER]. 
3. Program Options 
After the PCRSIM command has been issued to start 


the simulation, the first screen seen by the operator is at 


QZ 


Figure 55. The second screen the operator will see is at 


Figure 56. At this point the operator must choose which 


——————— 








This program will demonstrate the effect of various 
waveforms on the pulse compression process. The presence 
of a doppler shift will have an effect on the pulse 
compression process. The waveforms are implemented in one of 
three ways. The first way is to view the pulse compression | 
process with no doppler effects present. This will give the 
theoretical results of a best case scenario. The second way 
'is to input the effects of doppler by entering either the 
_doppler frequency or the maximum allowable phase shift you 
wish the waveform to undergo. The third method is to enter 
the effects of doppler and then use one of the standard 
correction schemes to try and cancel its effects. 











PRESS ANY KEY TO CONTINUE 







Figure 55: First Screen of PCRSIM main control program 


type of doppler processing the modulation scheme is to 
undergo. 

Once this has been done the main program transfers 
control to the subprogram selected in Figure 56. All of the 
screens for the subprograms are very similar, the only 
difference being what are the required inputs. 

The screen at Figure 57 shows the selection choice 
for the available modulation schemes if the no doppler 
choice was selected in Figure 56. The additional inputs 
required for the doppler or corrected doppler subprograms 
will be discussed further on. Once one of the modulation 


schemes has been selected than the required inputs will be 
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Which doppler technique do you wish to use? : 









To see pulse compression without doppler, ENTER 1 





To see pulse compression with doppler, ENTER @ 
To see pulse compression with corrected doppler, ENTER 3 
To exit the simulation, ENTER 0 


PRESS ANY KEY 20 COnNi ive 


Figure 56: Second Screen of PCRSIM main control program 


You have chosen not to consider doppler sensitivity. 

You now need to choose the waveform you wish to process. 
For a SINGLE BARKER CODE, choose "n=1", or 
For a COMBINED "m in n" BARKER CODE, choose "n=2Ziem 
For a Pseudo Random code, choose "n=3", or 
For a Polyphase code, choose "n=4", or 

For a Frank code, choose! n—-5 sc 

To return to the main menu, chnoese "n=6" > on 

To quit, choose "n=0" 


PRESS ANY KEY TO CONTINUE 


ENTER the # of your waveform choice, n = 





Figure 57 : The no doppler modulatiom ssencmogeeu--7 
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Peqguested for thatecode as discussed in Section A of this 


chapter. The screen for a Barker code is at Figure 58. 











— = 





The known Barker codes have lengths of n= 2,3,4,5,7,11 or 
13. You need to choose the length of the Single Barker code 
you want. 






Then press RETURN 






ENTER the # of the Single Barker Code to use, 






Figure 58 : The single Barker modulation scheme screen 

The results will be as for Figure 48 if a single 
Barker code of length 13 is selected. The next option is to 
use the eclipsing loss subprogram from the choices presented 


at Figure 59. 


If you want to do a Single Barker Code again, choose n 









Memvou Want to return to the non-doppler menu, choose n = 2 









If you want to return to the Main menu, choose n = 3 





Tf you want to view the PRI and pulse length, choose n = 4 











me yOu want to quit, choose n = 0 





moien press RETURN 











Choose n = 


Figure 59 : Final screen of a modulation scheme subprogram 


If the ECLIPS.m subprogram has been chosen then the 
screen at Figure 26 will request the necessary inputs. An 
example of the output of this subprogram is at Figure 27 for 
meee, PRE = 5 KHZ and 7 = 5 cycles of the RF for a PCR 


= 169. This plot shows the pulse width as compared to the 


iOS 


You will need to input values for the RF frequency 
of operation of the system, the PRF you wish to consider 
and the # of cycles of the RF that one Tau is to be 
comprised of. The PCR will either be given by a higher 
level program Ob (inst by eyecur 


Fnter the REF freq of Operation in 9) Mhz eee ee 
Enter the PRF in [ KHz }], PRF = 

Enter the # of cycles of RF that is 1 Tau, # = 
Enter 1 if PCR present or 2 £6 1INpUESIE, "eloteen— 


Figure 60 : First screen of the ECLIPS.m subprogram 


SP4CE PCR TAKES UP BETWEEN PRI, TAU tokes up 3S cycles of RFE 








= | 

= 1.3 7 

= 

Gye 
© 4 =z 3 + = 
PRI Imtervals for 3 pulses. in € sec J x ian 

Figure 61 : Eclipsing loss graphical presentation over 3 
PRI 


IONS 


PRI for three pulses. The unambiguous range and the minimum 
range are calculated and displayed for these inputs. After 
this the operator has the option to return to the main menu 
@eamieo redo the eclipsing loss option. 

Should any data entry error have occurred that the 
inherent MATLAB error routines did not recognize, then the 
Simulation will detect it and the error message at Figure 62 
will be displayed. All variables will be cleared, and the 


main program will be restarted. 


YOUR CHOICE DOES NOT CORRESPOND TO THOSE AVAILABLE. 


THE PROGRAM WILL BE RESTARTED FOR YOU. 


MeMeMpeR YOU MUST CHOOSE O TO EXIT. 


Beeoo ENTER TO CONTINUE... 





Figure 62: Error message subprogram screen 

The additional input screens necessary for the doppler 
pulse compression options are at Figures 63, 64 and 65. 
These screens appear after the modulation choice has been 
made because they need to use some of the inputs from the 
modulation scheme. The first screen of the doppler 
sensitivity choices available is at Figure 63. The choice 
is between entering the doppler frequency for a specific 
case or entering the overall phase shift to test a code's 
response to any system that could have a doppler velocity 


which would produce such a phase shift. 


ne 7 


The additional input screens necessary for the 


corrected doppler pulse compression options are at Figuimes 











You must now decide whether you want to see the 
result for a doppler frequency, or for a Maximum pies 
shift across the coding sequence. 








By giving the maximum phase shift, you will be 
controlling the amount of error in the processed signal. 
With this option, your maximum doppler will be dictated. 






By giving a doppler frequency, 
caused by that doppler. 


| Enter 1 to specify doppler shift, or 2 for phase Grmgume | 


Figure 63 : Doppler choice screen of the DOPP.m subprogram 


you will see the error 





Enter dopplemn freq. sn nHz 


Enter the compressed pulse width tau in microsec. 





Figure 64 : Doppler choice screen if doppler shift selected 


Enter the maximum allowable phase shift in radians 


Figure 65 : Doppler choice screen if phase error selected 





66 and 67. For the first screen at Figure 66, the pOimiiee 


note is that a positive doppler velocity entry means that 
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the correct decision as to whether it is an opening or 
closing target was made. A negative entry means that an 
opening target was identified as being a closing target or 


vice versa, and thus the consequences of an error will be 


Soccm. 





You have chosen to use the Doppler sensitivity 
feature. The simulation will assume that the correct 
decision as to whether it is an opening or closing target 
has been made if the doppler frequency entered is 

positive. If the doppler frequency entered is negative 
then you will see the consequences of having made 
the wrong decision. 









You must now enter the doppler data the program will 
require. | 






Figure 66 : Doppler choice screen if XDOPP.m selected 


Enter target maximum speed in knots 

Enter actual target speed in knots 

Enter the RF operating frequency in MHZ 

Enter the PRF in KHZ 

Enter the compressed pulse width in microseconds 


Compensation for the doppler shift by "derotating" 
using an assumed derotation doppler of one-half target 
maximum speed will be done. 





Figure 67 : Data entry screen if XDOPP.m selected 


E. INTERPRETATION OF RESULTS 

The output of the simulation is mostly graphical in 
nature Since graphical displays are readily assimilated and 
provide a more intuitive understanding of the process. 


Results from each of the subprograms will be presented. 
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These plots will be analyzed to demonstrate some of the 
Capabilities of the simulation program and to provide worked 
examples as a reference guide. 
1. Pulse Compression Without Doppler 

The subprograms without doppler are all quite self 
explanatory, and examples for each modulation scheme are at 
Figures 48, 49, 50, 51 and 52. This selection allows ume 
operator to select a coding scheme that meets his ideal 
needs, and the other two options allow him to see how this 
code performs under doppler conditions and if a correction 
scheme can help its performance. 

2. Pulse Compression With Doppler 

The subprograms with doppler allow that either a 
specific case (doppler frequency specified) or a general 
case (maximum allowable phase shift specified) be 
investigated. A combined 13 in 13 Barker that has undergone 
a maximum phase shift of 90° is at Figure 68. The output of 
the correlator has degenerated so much as compared to the 
ideal case at Figure 49 that it is useless for detection 
purposes. A P3 polyphase scheme with the same PCR as that 
of the combined Barker is at Figure 69. Here the output of 
the correlator is very Similar to that of the ideal case at 
Figure 51. Note however that a range error has been 


introduced as the mainlobe of the response of the correlator 


tgs) 


Combined 13 in 13 Borkxer Code wY doeppier 





Mox Prose SHiiftt = 1.571 Rods 
1S0O fF 


|Rotio of Doppler pulse 
! 
Somnipre ssiom Ca MCR a 3Si9 13 = 


SoG | 


Magnitude of the CROSS CORRELATION 





— 1 Ste, = fee; — 32 © S0° 100 1S So 


PCR = 169. PSL @» —4.624¢ 098 ond ISL = 3.786 coh 


Figure 68 : 13 in 13 Combined Barker with 90° of phase shift 
overall 


Psiyepnase Code with doppier for n = 169 
Soo Mox PrRose Shift = 1.971 Rods a 
| 
= 120 FRot!lo cof Ocppier puise 4 
— | 
| 
Pe = “y 
= 120 ormmpresslon to CR 90.03 *% 7 
a= 
—7 
cama 
wn 100 
a 
=) 
x= 
a> 
= 
ara] 60 
| 
= 
ce 420 
= 
[=] 
= 
FAS, 
Qa LE SE SSE EEE ete ee a nn ae Manet, Mala 


—150 —19090 —So OQ so 100 1S0O 


PCR = 169, PSL = —9.423 5S8 ond ISL = —5.732 9B 


Figure 69 : 169 component P3 polyphase code with 90° of phase 
shift overall 
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has been shifted in time. If the doppler shift that caused 
this is known, then the range error can be corrected. 
3. Pulse Compression With Corrected Doppler 
The subprograms with corrected doppler are always 

done for a specific case. A combined 13 in 13 Barker that 
has undergone a doppler phase shift is used in the 
correction scheme with the following parameters: 

® RF=2 GHZ 

@ PRE—= 130hnzZ 

® Tau = 1 psec 


@ V,, = 700 KPH 


tet 
® Vv. = 900 KPH. 
The output is at Figure 70, and clearly this simple doppler 
correction scheme was very effective. The output at Figure 
71 shows what the result would have been without the doppler 
correction being applied. This was done simply by letting 
Vioxx = O which gives the report for uncorrected doppler. The 
response at Figure 70 is a useable compression, albeit 
somewhat less than ideal. The correlator output at Figure 
71 would likely result in two targets at different ranges 
being declared rather than the one target that is actually 
there. The correlator output at Figure 72 is what would 
result if a wrong decision as to an opening or closing 


target had been made. This was done by entering £7) >see 


Gsael 


sey 


Combined 13 in 13 Borker Code with corrected doppler effects 


So FPatio of Ooeppler pulse | 


124Orcompression to PCR = 8452 &% 4 
ro | 4 
1004 = 





BOP 





Magnitude of the CROSS CORRELATION 





— 1» = Vere; 2; Oo He; 1090 US, 


Tirmme Cim Tau), PCR = 169, PSui =» —-19.77 3dB end iSul = -—-7 797 2B 


Figure 70 : Corrected doppler after pulse compression for a 13 
in 13 combined Barker code. The pulse compression achieved is 
Soe OL the PCR. 


Combined 13 in 13 Borker Code with corrected doprvrler effects 

169 FRatia of Doppler pulse a 
= 140 ormpreasaalom to PCR — 4144 % | 
=> | 
= 
o 120 
a= 
oO 
eS ole) 
wo 
Ss 
4 So 7 
a} 
w 
= 60 a 
s | 
3 20 - 
= | 
= 
S 20 - 
=z | 


—1350 —100 —So oO So 100 TSO 


Tirme Cim Tou), PCR = 169, PSl — —-0.3478 98 ond ISL @= 2458 3B 


Figure 71 : Uncorrected doppler after pulse compression fora 
13 in 13 combined Barker code. The pulse compression achieved 
1s 41% of the PCR. 


a bes 


C§ombpined 13 im 13 Barker Cade with corrected doppler effects 


Ratian of Ooppler pulse 
i =e) 


compression to PCR = 43.83% 


Magnitude of the CROSS CORRELATION 


—-150 —- 100 — =e, Oo (2) 1) OXe, eer 


Time Cim Tou), PCR = 169, PSL —2 611 98 ond ISL = 2634 2B 


Figure 72 : Corrected doppler after pulse compressionwiomae 
13 in 13 combined Barker code where an incorrect decision as 
to whether it 1S an opening or closing target was made. The 
pulse compression achieved is 43% of the PCR. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

This thesis has shown the background to the development 
of the DFT as a filter bank for radar applications. This 
development was then simulated in MATLAB for single filters, 
filter banks and the generation of target reports. In order 
to reduce the sidelobe response of the filters, weighted 
inputs through the use of windows were developed and 
Simulated for a weighted DFT filter bank. Finally, the 
concept of pulse compression was addressed and various 
representative pulse compression schemes were simulated 
under a variety of doppler conditions. 

Traditional classroom methods are not sufficient to 
demonstrate radar signal processing beyond trivial cases. 
The use of Simulation allows the student the freedom to go 
beyond the trivial cases and arrive at a fuller 
understanding of the principles involved. Simulation on a 
personal computer is simple, cost effective and accessible. 

These Simulation programs are freely available for 
student use and with them the computational burden of 
addressing some of the real world issues of radar signal 
processing has been removed. The simulations deliver the 


Same degree of accuracy as the real RSP processes would, 


Ves 


albeit at a much slower pace than the on board computerauas 
a modern radar. The simulation programs are available from: 
Professor G.S. Gill, 
Code EC/G1 
Naval Postgraduate school, 


Monterey, CA 93943 


B. RECOMMENDATIONS 


Some of the possibilities available for further thesis 
study in this area can be seen from the algorithm processing 
blocks in Figure 2. In particular, the simulaticonj@paeae 
processing blocks necessary for the adaptive nulling, CFAR 
detector and ambiguity resolution processes are considered 


to be a logical extension to this thesis. 
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